REMARKS 

Pursuant to 37 C.F.R. §1.111, reconsideration of the instant application, as amended 
herewith, is respectfully requested. Entry of the amendment is requested. 

Claims 1-1 1 1 are presently pending before the Office. Claims 12, 26, 51, 63, 65, 90, 92, 
93, 94 and 99 have been canceled, and claims 102-1 1 1 have been added. Applicant has amended 
claims by the addition of claims. No new matter has been added. Support for the amendments 
can be found throughout the specification as originally filed. Applicant is not intending in any 
mariner to narrow the scope of the originally filed claims. 

The Examiner's Action mailed November 10, 2003 (Paper No. 8) and the references cited 
therein have been carefully studied by Applicant and the undersigned counsel. The amendments 
appearing herein and these explanatory remarks are believed to be fully responsive to the Action. 
Accordingly, this important patent application is believed to be in condition for allowance. 

Applicant acknowledges and appreciates the request made by the Examiner to provide 
information as delineated in paragraph 2 of the Office Action under 37 CFR 1.105. Applicant 
herein submits his good faith Rule 132 affidavit attached herein, providing to the best of his 
knowledge, the information desired by the Examiner. 

Relying on 35 U.S.C. §1 12, second paragraph, the Office has rejected the subject matter 
of claims 1-101 as being indefinite for failing to particularly point out and distinctly claim the 
subject matter which Applicant regards as the invention. The Examiner alleges that it is unclear 
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what molecules, dimers and atoms are being claimed. Applicant respectfully traverses the 
rejection and requests reconsideration. 

Applicant submits that ALL kind of molecules, dimers and atoms can indeed form a 
magnecule without any restriction. 

Also, magnecules have been proved to exist at the gaseous, liquid and solid state. 
Alternatively, any gas magnecule can be brought to the solid state via cryogenic means during 
which the magnecular structure IS NOT lost because we decrease the temperature (not so for a 
magnecule at the solid state that can be lost under sufficient heat). 

Therefore, Applicant submits that claims 1-101 do define the legal metes and bounds of 
the invention. It is not the role of the claims to enable one skilled in the art to reproduce the 
invention but rather to define, for those skilled in the art the legal metes and bounds of the 
invention. 

It is respectfully submitted that claims 1-101 and new claims 102-121 fully comply with 
35 U.S.C. §112, second paragraph. Withdrawal of the rejection is respectfully requested. 

Relying on 35 U.S.C. §103(a), the Examiner has rejected the subject matter of claims 1- 
15, 18-21, 29-45, 52-67, 79-95 and 101 as obvious over Richardson (the '656 patent). Applicant 
respectfully traverses the rejection and requests reconsideration. 

Applicant herein incorporates by reference the attached Affidavit of Applicant Ruggero 
Maria Santilli submitted under 37 CFR 1.131 and 1.132, wherein applicant submits that the 
parent application predates the filing of the Richardson patent, and that Richardson admits that 
an eminent scientist is the inventor of the magnecule and implicitly states that it was discovered 
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before the filing of his application. Further, Applicant provides an affidavit of Leon Toups of 
Toups Technology wherein he verifies that Applicant was in possession of the invention prior to 
Mr. Richardson and that Mr. Richardson was in an employee position that allowed him access to 
the evidence developed by Applicant. 

Accordingly, the Examiner has not established a prima facie case of obviousness. 

Applicant notes that pages 15 and 34 of the present application refers to Articles co- 
written by Applicant and Donald D. Shillady of Virginia Commonwealth University, Department 
of Chemistry, Richmond, Virginia. Although the specification itself discloses the relevant 
portion, the complete arrticles are attached herein for the Examiner's perusal and consideration 



Applicant respectfully submits that claims 1-111 are patentable over the art of record. 
A Notice of Allowance is earnestly solicited. 

If the Office is not fully persuaded as to the merits of Applicant's position, or if an 
Examiner's Amendment would place the pending claims in condition for allowance, a telephone 
call to the undersigned at (727) 538-3800 would be appreciated. 



CONCLUSION 



Very respectfully, 





Dennis G. LaPointe 
Mason Law, P. A. 
17757 U.S. Hwy. 19 R, Suite 500 
Clearwater, FL 33764 



(727) 538-3800 
Reg. No. 40,693 
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Abstract 

Despite outstanding advances throughout this century, we still lack final knowledge on the structure 'of the hydrogen 
molecule because of a number of insufficiencies in available models identified in the text. In this paper we use the recently 
achieved covering of quantum chemistry called hadronic chemistry in its iso-chemical branch, and introduce a new 
model of the hydrogen molecule characterized by a bond at short distances of the two valence electrons into a singlet 
quasi-particle state called iso-electronium. We study the iso-chemical model of the hydrogen molecule with a stable iso- 
electronium describing an oo-shaped orbit around the respective two nuclei, and another model with a weaker realization 
of the iso-electronium as a partially stable state. We show that the new model provides, apparently for the first time, an 
exact representation of the binding energy and other characteristics of the hydrogen molecules from axiomatic principles, 
w ithout ad hoc modifications of theory. In subsequent papers we shall show that the new model permits apparently 
novel advances in the energy, liquefaction and other technological applications of the hydrogen. © 1999 International 
Association for Hydrogen Energy. Published by Elsevier Science Ltd. All rights reserved. 



1. Limitations of quantum chemistry 

Unquestionably, quantum chemistry [1-4] constitutes 
one of the most important scientific achievements of this 
century whose multiple and diversified applications have 
produced clear benefits to all of mankind. 

Nevertheless, science will never admit 'final theories', 
and structural generalizations of pre-existing doctrines 
are only a matter of time. 

Doubts on the final character of quantum mechanics, 
and, therefore, of quantum chemistry, can be traced back 
to Einstein. Podolsky and Rosen (E-P-R) [5]. who intro- 
duced their celebrated argument of 'lack of completion' 
of the theory. 

With the advancement of experimental and tech- 
nological knowledge, more and more limitations of quan- 
tum chemistry have been identified, such as 

( 1 ) Quantum chemistry does not characterize an attract- 
ive force among the neutral atoms of the H- (and 
other) molecule(s) sufficient to explain their bond. In 
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fact, the total Coulomb bonding force between atoms 
is identically null at the semi-classical level. For this 
reason, the bond is assumed to be due to exchange, 
van der Waals and other forces originating in nuclear 
physics whose strength is known to be weak. To 
state it differently, quantum chemistry still lacks in a 
molecular bond the equivalent of the strong force in 
nuclear physics. 

(2) Current models of molecular bonds do not explain 
why the hydrogen (and water) molecules admit only 
two H-atoms, and not three or more. In fact, the 
exchange, van der Waals and other forces were con- 
structed for the specific purpose of admitting an arbi- 
trary number of constituents as necessary to represent 
nuclei. This main characteristic evidently carries over 
at the level of molecular bonds, thus preventing a 
scientific explanation why nature has preferred two 
hydrogen atoms in the hydrogen and water molecule. 

(3) A number of main characteristics of H- (and other) 
molecules have not been represented accurately so 
far by quantum chemistry. This is the case for the 
binding energy, electric and magnetic dipole and 
multi-pole moments, liquid state and other aspects 
which still miss a few percentages in their represen- 
tation. 
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(4) Quantum chemistry does not allow meaningful ther- 
modynamical calculations. As an example, the 2% 
currently missing in the representation of molecular 
binding energies corresponds to about 959 kcal/mole, 
which is too large an uncertainty for chemical ther- 
modynamics where reaction energies are typically of 
the order of 20 kcal/mole. 

(5) More accurate numerical representations of molec- 
ular data can be achieved via the so-called 'screening 
of the Coulomb potential* which however constitute 
a structural departure of the theory from the basic 
axioms of quantum mechanics and chemistry. As an 
example, Gaussian screening of the Coulomb law 
generally requires non-unitary transforms of the orig- 
inal Coulomb settings, thus implying the exiting from 
the class of equivalence of quantum mechanics and 
chemistry. 

(6) Current models of molecular bonds imply a behavior 
of H- (and other) molecules under external electric 
and magnetic fields contrary to experimental, 
evidence. As an illustration, the atoms of a molecular 
bond preserve their individuality according to quan- 
tum chemistry. One of the best established disciplines 
of this century, quantum electrodynamics, then 
implies that the valence electrons of all atoms of a 
molecule must assume the same polarization under 
an external magnetic field. In turn, this implies the 
capability by water to acquire a paramagnetic charac- 
ter which is contrary to evidence. 

(7) Despite the availability of more powerful computers, 
accurate configuration interactions converge slowly, 
thus requiring excessive amounts of computer time 
for calculations. Such amounts of time increase with 
accuracy, reaching prohibitive amounts in attempt- 
ing exact representations. 

(8) The still unexplained 'correlation energy* is at times 
advocated to treat the missing percentages; orbital 
theories work well at qualitative and semi-empirical 
levels, but they remain afflicted by yet unresolved 
problems, such as the correlation among many elec- 
trons as compared to the evidence that the correlation 
occurs mainly for electron pairs. 

(9) When models achieve 100% representation of the 
experimental data, this occurs via the introduction of 
a number of empirical parameters which generally 
lack physical meaning and motivation; and other 
problems. 

All the above limitations refer to molecular structures 
which are closed-isolated and reversible (i.e., they verify 
usual total conservation laws and their time reversal 
image is causal). 

The limitations of quantum chemistry are multiplied 
when considering irreversible structures, such as chemical 
reactions involving hydrogen. In fact, we have in this case 
a structural incompatibility between the basic axioms of 
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quantum chemistry, and the reality represented because 
said axioms are structurally reversible (as a necessary 
condition to represent molecules), without any known 
possibility to modify them in such a way as to reach a 
credible representation of the irreversibility in chemical 
reactions from first axiomatic principles. 

When passing to more complex biological structures, 
any claim of final character of quantum mechanics and 
chemistry implies the exiling of the boundaries of science. 
In fact, quantum mechanics and chemistry can only rep- 
resent systems that are closed and isolated from the rest 
of the universe (conservative), reversible and rigid, as well 
known. On the contrary, biological systems require the 
representation of the time rate of variation of charac- 
teristics (e.g., their increase at birth, and their decrease at 
a later age), they are irreversible and flexible. To put it 
bluntly, under the conditions of the exact validity of 
quantum mechanics and chemistry everywhere for the 
microcosm, our body should be perfectly isolated from 
the rest of the universe, perfectly eternal and perfectly 
rigid (for studies of these biological aspects one may 
consult Ref. [15]). 

Moreover, biological systems are generally multi- 
valued, in the sense of requiring a structurally novel 
mathematics expectedly of the so-called hyper-structural 
type. In any case, the belief of understanding the 
extremely complex code of a DNA via the single-valued 
mathematics underlying quantum chemistry is beyond 
scientific credibility. 

In view of these and other occurrences, in this paper 
we introduce a covering of quantum chemistry which 
admits the latter as a particular case, yet includes new 
properties capable of resolving at least some of the above 
insufficiencies. 

We shall then introduce a new model for the structure 
of the hydrogen molecule which, apparently for the first 
lime, can: explain why the molecule has only two atoms; 
achieve an exact representation of the binding energy 
and other characteristics; permit a dramatic reduction of 
computer time; and reach other advances. 

In subsequent papers we plan to indicate that, even 
though admittedly small, the percentages currently miss- 
ing in the exact representation of the characteristics of the 
hydrogen molecule imply rather important technological 
advances, including possible new structures of the gas 
itself, new enhancement of the energy content, new means 
for liquefaction and other applications. 

2. Rudiments of hadronic mechanics 

In this paper we follow the historical E-P-R teaching 
on the Mack of completion* of quantum mechanics and, 
therefore, of quantum chemistry [5-9]. 

Among a variety of possible "completions*, we study 
the addition of effects resulting from the deep overlapping 
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of the wavepackets of valence electrons in chemical 
bonds. 

As it is well known, all massive particles (beginning 
with the valence electrons) have a wavepacket of the 
order of 1 fm (10~ 13 cm). We therefore assume quantum 
chemistry as being exactly valid at distances bigger than 
1 fm, while we seek possible, generally small corrections 
at mutual distances of 1 fm or less due to the deep over- 
lappings of valence electrons. 

Deep mutual penetrations of the wavepackets of par- 
ticles are known to be: non-linear (i.e., dependent on 
powers of the wavefunctions greater than one), non-local- 
integral (i.e., dependent on integrals over the volume of 
overlapping which, as such, cannot be reduced to a finite 
set of isolated points), and non-potential (i.e., consisting 
of 'contact* interactions with consequential 'zero range' 
for which the notion of action-at-a-distance. potential 
energy has no mathematical or physical sense). 

The latter features evidently imply a non-Hamiltonian 
theory (i.e., a theory which cannot be solely characterized 
by the Hamiltonian thus requiring additional terms). It 
then follows that the emerging theory is non-unitary (i.e., 
its time evolution violates the unitarity conditions 
Ulf = lfl) = I when formulated on conventional Hil- 
bert spaces over conventional fields). 

In short, we study a 'completion' of quantum mech- 
anics and chemistry via the addition of effects at distances 
of the order of 1 fm (only) which are assumed to be non- 
linear, non-local, non-potential, and non-unitary, here 
generically referred to as 'non-Hamiltonian effects*. Note 
that the condition of non-unitarity is necessary to exit 
the class of equivalence of quantum chemistry. 

As it is well known by experts in the field, the main 
technical difficulty of the above completion is that, when 
treated via the conventional formalism of quantum chem- 
istry, said short range non-Hamiltcnian effects have 
physical inconsistencies so serious to prevent any con- 
sistent application to the real world [11,16]. 

To give an idea, we indicate that: non-linear theories 
violate the superposition principle, thus having no mean- 
ingful application to composite systems such as mol- 
ecules; non-local theories violate causality; non-potential 
theories violate probability laws; non-Hamiltonian the- 
ories do not admit invariant basic units of time, length, 
energy, etc., thus having no meaningful applications to 
real measurements; non-unitary theories do not preserve 
the original hermiticity of operators, thus having no 
physically acceptable observables; and other incon- 
sistencies (see Refs [1 1, 16] for details). 

Quantitative studies on a non-Hamiltonian completion 
of quantum mechanics, capable of resolving the above 
problematic aspects, have been recently permitted by 
hadronic mechanics [10-21], which is an image of quan- 
tum mechanics characterized by the novel iso-, geno- and 
hyper-mathematics [10] for the representation of revers- 
ible, irreversible and multi-valued systems, respectively. 
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with Hamiltonian as well as non-Hamiltonian effects. It 
then follows that hadronic mechanics has three cor- 
responding branches called iso-mechanics, geno-mech- 
anics and hyper-mechanics, respectively. 

The latter new mathematics are essentially charac- 
terized by the lifting of the conventional unit / = + 1 
of quantum mechanics into invertible generalized units 
characterized by an n x n matrix (or operator) / with an 
arbitrary functional dependence, / = + 1 -+J — I 
{i.r.t .ilt.Ptp ,...) =J/T, which are Hermitean. /= /. 
non-Hermitean, /V /. and multi-vaiued-non-Hermitean, 
/ = (/,. /, respectively. 

The new mathematics then emerge via the recon- 
struction of all conventional algorithms with respect to 
said generalized units /. thus including corresponding 
lifting of numbers and fields, metric and Hilbert spaces, 
algebras and geometries, etc. For instance, all con- 
ventional (associative) products A x B of quantum mech- 
anics have to be lifted into the generalized products 
A k B = A x fx Sunder which / = Iff is indeed the cor- 
rect left and right new unit, f*A = Axl'=A. called iso-, 
geno- and hyper-units. Similar generalizations hold for 
all other operations, including differentials and deriva- 
tives, conventional and special transforms, etc. [15]. 

The mam physical reason for the lifting of the trivial 
unit 7^+1 into an operator / is that it provides the 
only known invariant representation of non-Hamiltonian 
effects in the pairing of valence electrons [16J. 

To state it differently, in order to prevent possible 
misrepresentations of reality, contact-non-potential 
effects should be represented with anything except a 
(Hermitean) Hamiltonian. All representations of such 
effects other than that via Hamiltonians have been proved 
to be non-variant, thus unsuitable for applications. The 
representation of contact-non-potential effects via gen- 
eralized units has been selected because it is the sole 
known approach yielding an invariant theory which, as 
such, admits predictions suitable for experimental veri- 
fications. 

As a simple example, the representation of the novel 
interactions caused by deep wave-overlappings of the 
valence electrons as in Fig. 1 can be only done in quantum 
mechanics and chemistry via the addition of potentials 
in the Hamiltonian. However, this implies granting 
potential energy to contact effects which have none (it 
would be like representing the drag force experienced by 
the space-shuttle during re-entry with a Hamiltonian), 
thus implying the description of a system that has no 
connection with physical reality. 

A consistent representation of the same non-Ham- 
iltonian interaction of Fig. 1 can be easily achieved by 
hadronic mechanics via the generalized unit called iso- 
unit 

/=cxp[Cfd-V^(r«f 1 (r)], (2.1) 
where C is a constant (at short range) and ^ are the 
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Fig. 1. A schematic view of the fundamental conditions studied 
in this paper, the deep overlapping of the wavepackets of valence 
electrons in singlet coupling, as required by Pauli's principle, 
whose non-linear, non-local and non-potential structure is 
beyond the capabilities of quantitative treatment via quantum 
mechanics and chemistry, thus requesting a covering theory. 



wavef unctions of the two valence electrons evidently in 
singlet couplings (from Pauli's principle). Achievement 
of invariance then requires the step-by-step recon- 
struction of the entire formalism of quantum mechanics 
with respect to the new unit (2.1). 

Hadronic representation (2.1) of the new interactions 
of Fig. 1 has a number of advantages over the quantum 
form, such as: the treatment of non-potential effects with- 
out a potential; the achievement of an invariant for- 
mulation; the preservation of hermiticity-observability 
under non-Hamiltonian effects: and others [14, 15]. 

Whenever the wave-overlapping is no longer appreci- 
able, the exponent of generalized unit (2.1) is null. Then 
7=+l, and hadronic mechanics recovers quantum 
mechanics identically at all levels. This illustrates the 
covering nature of hadronic over quantum mechanics, 

A similar situation occurs for other characteristics 
which are conceptually, mathematically and physically 
beyond the representation capabilities of a Hamiltonian, 
but which admit instead a simple and effective rep- 
resentation via generalized units. 

An illustration is given by the representation of actual, 
extended, non-spherical and deformable shapes of bodies 
which is impossible for quantum mechanics because it 
would violate one of its pillars (the rotational symmetry), 
but which can be readily represented by hadronic mech- 
anics. As an example, a deformable spheroidal ellipsoid 
can be represented via the iso-unit /= Diag. n]), 
where the /f s are the variable semi-axes, and the recon- 
struction of the rotational symmetry with respect to that 
generalized unit (see Refs [9, 15] for details). 

Similarly, Bose-Einstein correlations and con- 
densations are known to be beyond the exact character 



of quantum mechanics because intrinsically extended in 
space, thus requiring a suitable non-local topology. Had- 
ronic mechanics can readily represent these effects from 
first axiomatic principles via annxn generalized unit 
incorporating off-diagonal elements of type (2.1). 

We refer the interested reader to Ref. [1 1] for an outline 
of numerous applications and experimental verifications 
of hadronic mechanics in particle physics, nuclear phys- 
ics, superconductivity, astrophysics and other fields. 

A mathematical motivation for the generalization of 
the unit is that it permits the reconstruction of the linear, 
local and potential characters on certain generalized 
spaces (called iso-, geno- and hyper-spaces) over gen- 
eralized fields (called iso-, geno- and hyper-fields). 

The latter features assure the resolution of the prob- 
lematic aspects indicated earlier for the conventional 
treatment of non-Hamiltonian effects. The same features 
also imply that hadronic mechanics coincides with quan- 
tum mechanics at the abstract, realization-free level in all 
its iso-, geno- and hyper-branches, and merely provides 
broader realizations of conventional quantum axioms. 

In particular, hadronic mechanics preserve all con- 
ventional laws and principles of quantum mechanics, 
such as Pauli's exclusion principle, Heisenberg's uncer- 
tainties, causality and probability laws, etc. 

The understanding of this note requires a knowledge of 
hadronic mechanics and its underlying new mathematics 
that we cannot possibly outline here. Readers are there- 
fore referred to Ref. [10] on the new mathematics, Ref. 
[1 1] on the latest formulation of hadronic mechanics, and 
Ref. [16] on the inconsistencies of conventional for- 
mulations of non-Hamiltonian effects. 



3. Rudiments of hadronic chemistry 

In this paper we use the covering of quantum chemistry 
characterized by hadronic mechanics, here called had- 
ronic chemistry. As such, the new discipline too has three 
main branches called iso-, geno- and hyper-chemistry for 
the representation of reversible, irreversible and multi- 
valued chemical structures, respectively, with Ham- 
iltonian (long range) and non-Hamiltonian (short range) 
effects. 

Since the hydrogen molecule is stable and reversible, 
in this paper we shall focus solely on iso-chemistry. 

However, the reader should be aware that, as it 
occurred for conventional chemistry, iso-chemistry too 
is insufficient to represent chemical reactions, or, more 
generally, irreversible complex biological structures. 
When studying these cases in future papers the use of the 
broader geno-chemistry and hyper-chemistry, respec- 
tively, will be necessary for consistency. 

Again, by conception and construction, hadronic and 
quantum chemistry coincide everywhere, except at dis- 
tances of 1 fm or less in which the former admits novel 
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contributions over the latter under the exact validity of 
all conventional quantum laws. As we shall see, the novel 
short range contributions are generally small, yet they 
permit the apparent resolution of at least some of the 
insufficiencies outlined in Section 1 . 

The completion of quantum into hadronic mechanics 
along the historical E-P-R teaching is studied in detail 
in Ref. [9]. The same results hold for the completion of 
quantum into hadronic chemistry. For this reason, the 
topic is ignored hereon for brevity. 

Quantum mechanics and chemistry have only one for- 
mulation, the conventional one [1-4]. On the contrary, 
hadronic theories have two different formulations [10- 
21], that on generalized spaces over generalized fields and 
its projection on conventional spaces over conventional 
fields [10-16]. The reader should be aware that, to avoid 
excessive complexities for non-initiated readers, in this 
paper we shall solely consider the projection of hadronic 
theories on conventional spaces over conventional fields, 
and leave the mathematically correct formulation to 
mathematical studies. 

Moreover, systems are characterized by quantum 
mechanics and chemistry via the sole assignment of the 
Hamiltonian and the tacit assumption of the simplest 
possible unit / = + 1. By comparison, systems are char- 
acterized by hadronic mechanics and chemistry via the 
assignment of two quantities, the Hamiltonian for the 
representation of all conventional interactions, plus the 
generalized unit for the representation of the short range 
non-Hamiltonian effects. 

Despite apparent mathematical complexities, iso- 
chemistry can be uniquely and unambiguously con- 
structed in its entirety via a non-unitary transform 
characterizing precisely the iso-unit 



Uxlf = /V/, 



(3-D 



which must be applied to the totality of the conventional 
formalism as a necessary condition to achieve invariance 
[16]. 

This procedure yields the iso-unit itself, the iso- 
numbers, the iso-Hilbert space, etc. x ' 

j^UxIxlf = !=\/t f=(Uxlfy\ 
c-*Uxcxlf = cx(Uxlf) = cxf=c, 
AxB^Ux(AxB)xlf 

= (UxAxlf)x(Ux lf)~ y x(UxBxlf)~AxB, 

Hx\ii/y^ux(Hx\^y) 

= (UxHxlf)x{Uxlfy l x(Vx\\l,)) = Hx |f>, 
= «<fl x if) x ( V x Ify 1 x ( V x Wr» x(Uxlf) 
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(3.2) 



= <^|x|^> x/'ctc. 

Once achieved via the above simple map, iso-chemistry 
remains invariant under any other possible non-unitary 
transform, provided that it is evidently expressed in the 
iso-unitary form on the iso-Hilbert space over iso-fields, 
under which we have the invariance of the numerical 
value of the iso-unit, the invariance of the iso-product 
(where the numerical value of f is not changed), etc. 

H'xH* = /*/, W=Wxt\ 
Hx H* = W'xH* = H'x H>=/, 
/-/ = Wxix = /, 

A x tf- Wx {A x B) x l&* = A' xfxB' = A x B'.cXc. 

(3.3) 

The main problem of completing chemistry into iso- 
chemistry is, therefore, the identification of a transform 
which is: 

(1) non-unitary for r < 1 fm, Ux ifr * 1 fm = /> h 

(2) unitary for r > 1 fm, U x Ifr > 1 fm = 7; and 

(3) capable of representing the desired short-range, non- 
linear, non-local and non-potential effects due to the 
overlapping of the wavepackets of the valence elec- 
trons, e.g., as in iso-unit (2.1). 

All original axiomatic properties of quantum chemistry 
are preserved everywhere (including at short distances) 
and only realized in a broader way. This implies the 
preservation of all quantum laws and principles, such as 
Pauh's exclusion principle. Heisenberg uncertainty prin- 
ciple, causality and probability laws, etc. 

Regrettably, to remain within the limitation on length 
of this journal, we cannot provide additional details. A 
comprehensive presentation of hadronic chemistry, 
including its application to hydrogen, is forthcoming in 
Ref. [22]. 

We should finally mention that, by its very structure, 
iso-chemistry assures much faster convergence of power 
series. This can be illustrated by the fact that, any given 
divergent quantum series can always be turned into an 
isotopically convergent form via the selection of isotopic 
elements f sufficiently smaller than one [15], as in the 
example 



A(k) = A(Q) + k{A.Hll\!+- 



* > 1, 



A(k) = A(Q) + k[A 9 -H\!\.'+--- = N<co. f =/T\ 
[A,H] - AxB-BxA~{A r H] = A x H-HxA. (3.4) 

As we shall see, the above feature permits the reduction 
of computer time in chemical calculations by a factor of 
at least 1000, thus rendering feasible calculations that 
would be otherwise prohibitive. 
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Needless to say. the studies presented below are deeply 
linked to numerous preceding studies in quantum chem- 
istry (see, e.g., Refs [23-28]). In fact, we essentially pro- 
vide a 'completion 1 of available contributions via novel 
non-Hamiltonian effects at short distances without exit- 
ing f rom the axiomatic structure of the theory, as inherent 
in Coulomb screenings and other ad hoc modifications 
of quantum chemistry. We shall then show the achieve- 
ment of exact and invariant representations of main data 
of the hydrogen molecule. Other advances will be pre- 
sented in future papers. 



4. Iso-chemical molecular bonds 

We now present, apparently for the first time, the con- 
ceptual foundations of our iso-chemical model of molec- 
ular bonds for the simplest possible case of the H : 
molecule. We shall then extend the model to the water 
and to other molecules in subsequent papers. 

Since the nuclei of the two H-atoms remain at large 
mutual distances, the bond of the H : molecule is evidently 
due to the peripheral electrons, as generally acknowl- 
edged [1-4]. 

Our main assumption is that pairs of valence electrons 
from two different atoms bond themselves at short dis- 
tances into a singlet quasi-particle state, here called iso- 
electronium. which describes an oo-shaped orbit around 
the respective two nuclei, as il is the case of a planet in 
certain binary stars. 

We should stress that the above hypothesis is strictly 
prohibited by quantum mechanics and chemistry because 
of the repulsive Coulomb force which prevents any poss- 
ible bound state of identical electrons at short distances. 

However, as we shall see in the next section, the 
hypothesis is fully plausible for hadronic mechanics and 
chemistry because the novel non-Hamiltonian inter- 
actions due to deep wave-overlappings become so attract- 
ive w hen in singlet couplings at short distance, to over- 
come repulsive Coulomb forces. 

The above occurrence is important to see that the new 
model of the hydrogen molecule presented in this paper is 
strictly impossible for quantum mechanics and chemistry, 
and requires a covering theory. 

Iso-chemistry therefore introduces a new attractive 
force for the representation of the molecular bond among 
neutral atoms given by the coupling of valence electrons 
into singlet states. The binding energy is represented by 
the orbital motion of the iso-electronium around the 
respective nuclei. 

In this paper we show that the emerging new model of 
the hydrogen molecule can: provide an exact rep- 
resentation of the binding energy from first axiomatic 
principles: explain the reason why the hydrogen molecule 
has only two atoms: achieve a much faster convergence 
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of power series with consequential large reduction in 
computer times; and permit other advances. 

The explanation by iso-chemistry of the reasons why 
the H- (or water) molecule admit only two hydrogen 
atoms is inherent in the very notion of iso-electronium. 
In fact, once two valence electrons are bonded into the 
iso-electronium, there is no possibility for bonding 
additional atoms. 

The exact representation of the binding energy and 
other characteristics will be presented in the subsequent 
sections. 

Two notions of hadronic chemistry are important for 
further advances: 

(a) the hadronic horizon, which is the ideal sphere of 
radius 1 fm outside which quantum chemistry is 
assumed to be exact and within which hadronic chem- 
istry applies; and 

(b) the trigger, which is given by external (conventional) 
interactions causing identical electrons to move one 
toward the other and to penetrate the hadronic hor- 
izon against their repulsive Coulomb interactions. 
Once inside said horizon, the attractive hadronic for- 
ces overcome the repulsive Coulomb interaction, 
resulting in a bound state. 

The hypothesis of the bonding of electrons at short 
distances was first introduced by Sanlilli [29] in 1978 for 
the structure of the tt c meson as a hadronic bound state of 
one electron and one positron. Animalu [30] and Animalu 
and Santilli [31] extended the model to the Cooper pair 
in superconductivity as a hadronic bound state of two 
identical electrons. 

In the case of the tv model there is no need of a trigger 
because the constituents have opposite charges, thus nat- 
urally attracting each other. On the contrary, the existence 
of the Cooper pair does indeed require a trigger, which 
was identified by Animalu [30] and Animalu and Santilli 
[31] as being constituted by the Cuprate ions. For the case 
of an isolated hydrogen molecule, we conjecture that the 
trigger is constituted by the two H-nuclei which do indeed 
attract the electrons. We essentially argue that the attrac- 
tion of the electrons by the two nuclei is sufficient to cause 
the overlapping of the two wavepackets, thus triggering 
the electrons beyond the hadronic horizon. 

We should finally indicate that the hypothesis of the 
iso-electronium £an also be explained via the use of the 
intrinsic magnetic moments of the electrons which, in 
singlet couplings, are opposite to each other, thus imply- 
ing an attraction that can counterbalance the electric 
repulsion. 

However, the approach requires alteration (called 
mutation) of the electromagnetic fields of the electrons 
w r hich is directly representable by the Poincare-Sanlilli 
iso-symmetry [12, 13]; namely, an isomorphic recon- 
struction of the conventional Poincare symmetry with 
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respect to the iso-unit. This requires, for quantitative 
treatment, the lifting, this time, of relativistic quantum 
mechanics into relativistic hadronic mechanics which is 
here ignored for brevity [15]. 

The latter occurrence also explains that the term *iso- 
electroniunT is motivated by the fact that its constituents 
are not ordinary electrons (i.e., unitary representations of 
the Poincare symmetry), but rather electrons in mutated 
forms called iso-electrons (i.e., iso-unitary iso-rep- 
resentations of the Poincare-Santilli iso-symmetry). 



5. The stable iso-electronium 

We begin our quantitative analysis with the non-rela- 
tivistic quantum mechanical equation of two ordinary 
electrons in singlet couplings, e, and e T , represented by 
the wavefunction ip u {r) = \fr(r) and the conventional Sch- 
rodinger's equation 

H x ^(r) = x p + x </r(r) = E x 0(r). (5. 1 ) 

To transform this state into the iso-electronium rep- 
resenting the bonding of one valence electron with a 
valence electron of another atom of generic charge ze, we 
need first to submit eqn (5.1) to a non-unitary transform 
U x if = i # I characterizing the short range hadronic 
effects, and then we must add the trigger, namely; the 
Coulomb attraction by the nuclei. 

By recalling rules (3.2), this procedure yields the iso- 
Schrodinger equation for the iso-electronium (see [22] for 
details) 

UxHx\p(Pr) + trigger 

= UxHxlf x(Ux IT)- 1 x (/x^(r) + trigger 
= //x^(r)-h trigger 

= (-pxfpxf+-xf x <£(r) = £ 0 x ^(r), 

px\p{r) = -ixfxV^r), (5.2) 

where \p = t/xf H = UxHxlf, p= Uxpxlf, the 
factor f=(Uxlf)' [ in the first Coulomb term orig- 
inates from the non-unitary transform of model (5.1), 
while the same factor is absent in the second Coulomb 
term because the latter is long range, thus being con- 
ventional. 

The angular component of mode! (5.2) is conventional 
[1-4] and it is hereon ignored. For the radial component 
r = \r\ we assume the iso-unit [29-31] 

/= e w * * l+Afx^, . 

tf = JdrVV),*^ 

f * 1-TVxi^. (5.3) 
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Note that the explicit form of $ is of Coulomb type, 
thus behaving like Cxexp(-/>xr), with C approxi- 
mately constant at distances near the hadronic horizon 
of radius r r =l/£>, while ip behaves like 
D x ( 1 - exp ( — b x r))/r, with D being also approximately 
constant under the same range [29]. We then have 

e-"* r 

\-l A ir=\- K 0 . (5.4) 

(1-e~"")/> 

but the Hulten potential behaves at small distances as the 
Coulomb one: 

Vr* 1**^"x-. (5.5) 
b r 

As a result, inside the hadronic horizon we can ignore 
the repulsive Coulomb force altogether and write 

->->■► •» 
e- . :xe" e , :xc 

+ -xr- — a + - x(i - rV) 

/■ r r r 

e-* x ' 

= -K — , (5.6) 

by therefore resulting in the desired overall attractive 
force among the identical electrons inside the hadronic 
horizon, as desired. 

By assuming in the first approximation |fl = p * L 
the radial equation of model (5.2) reduces to the model 
of 7i meson or of the Cooper pair, although with different 
values of Fand b. which has been studied in great detail 
in Refs [29-31]. including all necessary boundary 
conditions. The solution is the typical one of the Hulten 
well: 

p-xfrxb 2 f mxV 1 V _ 

|£ u | = *— x - -n . (5.7) 

4xw {fxtfxb 2 » / 

To reach a numerical solution, we introduce the par- 
ametrization as in Ref. [29], k % = \//.xb, 
k 2 = m x Vjp 2 x ir x b 2 . We note again that, from the 
boundary conditions. k 2 must be bigger than but close to 
1 [29-31]. We therefore assume in first non-relativistic 
approximation that k 2 = m x Vjp 1 x h 2 x b 2 = 1 . By 
assuming that Kis of the order of magnitude of the total 
energy of the iso-electrons at rest as in the preceding 
models [29-3 1 ]. we reach the numerical value of the depth 
of the hadronic well (Hulten potential): 

|/* 2x/>xo>*2x0.5MeV = I MeV. (5.8) 

By recalling that p * 1, we reach the following numerical 
value of the radius of the hadronic horizon (radius of the 
iso-electronium): 

r t = /r 1 * (h 2 /m x F) 1 2 = (*/'» x «>„)' 2 
= (1.054 x 10" :7 ergsec 



• 
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.1.82xl0-- 7 g^><l■236xl0 2o Hz) , ' 2 

= 6.8432329 x 10" 11 cm = 0.01 5424288 bohrs 

= 0.006843 A. (5-9) 

It should be noted that: (1) the above value is only an 
upper boundary in the center-of-mass frame of the iso- 
electronium. i.e., it is the largest possible value under 
the assumptions of this section; (2) the value has been 
computed under the approximation of identically null 
relative kinetic energy of the iso-electrons with individual 
total energy equal to their rest energy; and (3) the value 
evidently decreases with the addition of the relative kin- 
etic energy of the iso-electrons (because this implies the 
increase of m in the denominator). 

The actual radius of the iso-electronium is also 
expected to vary with the trigger, that is, with the nuclear 
charges, as supported by the calculations presented in the 
next sections. This illustrates again the upper boundary 
character of value (5.9). 

The value k } is then given by A', = V/2xk 2 x 
bxc 0 = 0A9. Intriguingly, the combined two values 
fc, = 0.19 and k 2 = 1 for the iso-electronium are quite 
close to the corresponding values of the n° [29] and of 
the Cooper pair [30, 31]. 

It is important to see that, at this non-relativistic 
approximation, the binding energy of the iso-electronium 
is not only unique, but also identically null. 

The notion of a bound state with only one allowed 
energy level (called 'hadronic suppression of the atomic 
spectrum') is foreign to conventional quantum mechanics 
and chemistry, although it is of great importance for 
hadronic mechanics and chemistry. In fact, any excitation 
of the constituents, whether the n\ the Cooper pair or 
the iso-electronium ; causes their exiting the hadronic hor- 
izon, by, therefore re-acquiring the typical atomic spec- 
trum. Each of the considered three hadronic states has, 
therefore, only one possible energy level. 

In fact, it is easy to see that the generally finite Hulten 
spectrum reduces to only one energy value according to 

! * wi ~ 4xm \p2 xh 2 xb 2 ) 

= V x(ft ,_i)=«0. (5-10) 
4 x £ 2 

The additional notion of a bound state with null bind- 
ing energy is also foreign to quantum mechanics and 
chemistry, although it is another fundamental charac- 
teristic of hadronic mechanics and chemistry. In fact, the 
hadronic interactions admit no potential energy and, as 
such, they cannot admit any appreciable binding energy, 
as typical for ordinary contact zero-range forces of our 
macroscopic New tonian reality. 

The null value of the binding energy can be confirmed 
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from the expression of the meanlifeof the iso-electronium 
which can be written [29] 

T = fi/4x7rx/ 2 |^(0)|xa(xE^ n 

= 7.16 x 10 4 x ! = oo. (5.1 1) 

(k 2 — I) 3 *bxc 0 

The full stability of the iso-electronium, t = x, therefore, 
requires the exact value k 2 = 1 which, in turn, implies 
£ 0 = 0. 

According to the above non-relativistic model, the iso- 
electronium is a stable bound state (at ordinary 
conditions). This property is at the foundation of our 
representation of the stability of the H : and other 
molecules. We should, however, indicate that, under a 
relativistic treatment [11], the iso-electronium admits a 
small instability due to the exchange, van der Waals and 
other forces. 

By comparison, the Cooper pair is not permanently 
stable already at the non-relativistic level because its bind- 
ing energy is very small, yet finite [30], thus implying a 
large yet finite meanlife. Also by comparison, the n° can- 
not be stable, and actually has a very small meanlife, 
evidently because the constituents are a particle-anti- 
particle pair and, as such, they annihilate each other 
when bound at short distances. By comparison, no such 
annihilation may occur in the iso-electronium. 

Another important result of this section is that the iso- 
electronium is sufficiently small in size, as per value (5.9), 
to be treated as a single quasi-particle. This property 
will permit rather important simplifications in the iso- 
chemical structure of molecules studied in the following 
sections. 

By comparison, the Cooper pair has a size much bigger 
than that of the iso-electronium. This property is fun- 
damental in preventing the Cooper pair taking the role of 
the iso-electronium in molecular bonds, i.e., even though 
possessing the same constituents and similar physical ori- 
gins, the iso-electronium and the Cooper pair are differ- 
ent, non-interchangeable, hadronic bound states because 
they originate from different triggers and conditions. 

The lack of binding energy of the iso-electronium is 
perhaps the most important information of this section. 
In fact, it transfers the representation of the binding 
energy of molecular bonds to the motion of the iso- 
electronium in a molecular structure, as studied in the 
following sections. 

A novelty of iso-chemistry over quantum chemistry is 
that the mutual distance (charge diameter) between the 
two iso-electrons in the iso-electronium could, as a limit 
case, also be identically null, that is, the two iso-electrons 
are superimposed in a singlet state. Rather than being far 
fetched, this limit case is intriguing because it yields the 
value -2e for the charge of the iso-electronium, the null 
value of the relative kinetic energy, and an identically 
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null magnetic field. This is a perfectly diamagnelic state 
which evidently allows a better stability of the iso-chemi- 
cal bond as compared to a quasi-particle with non-null 
magnetic moment. 

Note that, if conventionally treated (i.e., represented 
on conventional spaces over conventional fields), the 
non-unitary image of model (5.1) would yield non- 
invariant numerical results which, as such, are unac- 
ceptable (Section 2 and Ref. [16]). This occurrence man- 
dates the use of the covering iso-chemistry and related 
iso-mathematics which assures the achievement of 
invariant results. 

Note also that the main physical idea of iso-unit (5.3) 
is the representation of the overlapping of the wavepack- 
ets of the two electrons under the condition of recovering 
conventional quantum chemistry identically, whenever 
such overlapping is no longer appreciable. In fact, for 
sufficiently large relative distances, the volume integral 
of iso-unit (5.3) is null, the exponential reduces to /, the 
non-unitary transform becomes conventionally unitary, 
and quantum chemistry is recovered identically. 

It is also important to see that, under transform (5.3), 
model (5.1) is implemented with interactions which are: 
non-linear, due to the factor ipfy in the exponent; non- 
local, because of the volume integral in (5.3); and non- 
potential, because it is not represented by a Hamiltonian. 

We finally note that the explicit form of the isotopic 
element f, eqn (5.3b), emerges in a rather natural (and 
unique) way as being smaller than one in absolute value, 
\f]« 1. This property alone is sufficient to guarantee 
that all slowly convergent series of quantum chemistry 
converge faster for iso-chemistry. 

6. Iso-chemical model of the hydrogen molecule with 
stable iso-electronium 

We are now sufficiently equipped to present, appar- 
ently for the first time, the iso-chemical model of the 
H : molecule. The understanding of the reasons why the 
molecule has only two H-atoms is inherent in the very 
concept of iso-electronium (Sections 4 and 5). In this 
section we shall, therefore, identify the equation of struc- 
ture of the H-molecule. The exact representation of the 
binding energy and other characteristics is studied in the 
subsequent sections. 

Our foundation is the conventional quantum model of 
the H : molecule [1-4]. The task is that of subjecting such 
a model to a transform which is non-unitary only at the 
short mutual distances r, = b~ 1 = r 12 of the two valence 
electrons (here assumed to be inside the hadronic 
horizon), and becomes unitary at bigger distances 
f r <; 10"" cm # /, I r » 10"" cm = /. This guarantees 
that our iso-chemical model coincides with the conven- 
tional model except for small contributions at r, = b~\ 

We assume that the state and related Hilbert space can 
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be factorized in the familiar form (in which each term is 
duly symmetrized or anti-symmetrized) as in Refs [1-4]. 

The non-unitary transform we are looking for shall act 
only on the r l2 variable, while leaving all others 
unchanged. The simplest possible solution is given by 
V(r i2 )xlf(r i2 ) = / = cxp([^(r 12 )/0(r l2 )]jdr 1 #(r l2 ) li x 
W'r)^), where the ^'s represent conventional wave- 
function and the ^s represent iso-wavefunctions. 

As a variant to the approach of Section 5 yielding the 
same end results, we construct the model by transforming 
short-range terms (iso-chemistry) and adding untrans- 
formed long-range ones (chemistry), thus resulting in the 
radial equation 

( - f x V, x f x V, - f x V. x f x V 2 + — 
V 2xji, ' 2x/i : " r l2 

xf_------- + ~W> = £x ^>- (6.1) 

*\a T 2a r 2h *<) 

By recalling that the Hulten potential behaves at small 
distances like the Coulomb one, eqn (6.1) then becomes 



- V 7 — \ 2 —v * r 

2x/i, 2x// 2 1— e" r «= x * 

_£l_^_^_-£l + ^W> = £x^>. (6.2) 
r Ul r 2a r lh r 2h R] 

The above equation does indeed achieve our objectives. 
In fact, it exhibits a new explicitly attractive force among 
the neutral atoms of the molecule which is absent in 
conventional quantum chemistry. The equation also 
explains the reasons why the H 2 molecule admits only 
two H-atoms (Section 5). As we shall see in the remaining 
sections, eqn (6.2) also permits the exact representation 
of the binding energy, yields much faster convergence of 
series with much reduced computer times, and resolves 
other insufficiencies of conventional models. 

Our iso-chemical model of the hydrogen molecule, 
eqns (6.2), can be subjected to an additional sim- 
plification which is impossible for quantum chemistry. 
In our isotopic model the two iso-electrons are bonded 
together into a single state we have called iso-electronium 
In particular, the charge radius of the latter is sufficiently 
small to permit that in first approximation, for r 12 « r Ut 
and r 1A , we have r Xa * r 2a = r a as well as, separately, 
r ]h * r 2h = r h . Moreover, the H -nuclei are about 2000 
times heavier than the iso-electronium. 

Therefore, our isotopic model of the H 2 molecule can 
be reduced to a restricted three-body problem similar to 
that possible for the conventional Ht ion [1-4], but not 
for the conventional H 2 molecule. It consists of two H- 
protons at rest and the iso-electronium moving around 
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Fig. 2. A schematic view of the proposed iso-chemical model of 
the hydrogen molecule as well as of other two-body molecules 
or. more generally, of the bond of one H-atom to a generic other 
element, e.g., HO, HC etc.. here assumed at absolute K and in 
the absence of rotations, vibrations, translations and other 
effects The fundamental assumption is that the two valence 
electrons, one per each atom, correlate themselves into a bonded 
singlet state at a short distance we have called *iso-electronium\ 
which is assumed in this case to be stable, w ith infinite lifetime 
as in eqn (5.9) and describing an oo-shaped orbit around the 
respective two nuclei. The iso-electronium is then responsible 
for the attractive force between the two atoms. The binding 
energy is instead characterized by the oo-shaped orbit of the iso- 
electronium around the two nuclei, conceptually represented 
in this figure via a standing wave. The model then permits a 
representation of the reason why the H : and H : 0 molecules have 
only two hydrogen atoms, the exact representation of the binding 
energy and other advances studied in the text. 

them in the oo-shaped orbit of Fig. 2, according to the 
structural equation 



tr fr . e r ' 

— Vt- — V;-Kx — — 
2fi i 2/n * j_ e - 



*^e~ *^e~ e~\ 
-=f -^ + ^W> = £x|^>. (6.3) 

Under the latter approximation, the model permits, 
for the first time, the achievement of an exact solution 
for the structure of the H-molecule, as it is the case for 
all restricted three-body problems. This exact solution 
will be studied elsewhere. In the next sections we shall 
study the solution of model (6.3) via conventional vari- 
ational methods. 

Note that the above exact solution of the hydrogen 
molecule is only possible for the case of the fully stable 
iso-electronium. In fact, for the case of the unstable iso- 
electronium the model is a full four-body structure which, 
as such, admits no exact solution. 

Note also that model (6.3) is the iso-chemical model of 
the H : molecule inside the hadronic horizon. The mat- 
ching representation outside the hadronic horizon is pre- 
sented below. Note also that the above restricted three- 
body model can be used for the study of the bonding of 
an H-atom to another generic atom, such as HO. thus 
permitting, again for the first time, novel exact cal- 
culations on the water as HOH: namely, as two inter- 
secting isotopic bonds HO and OH, with possible exten- 
sion to molecular chains and other molecules. 




X < Ifr" sec 



H + 



Fig. 3. A schematic view of the iso-chemical model of the hydro- 
gen molecule with unstable iso-electronium due to the weakening 
of the Hulten potential and consequential relaxation of the 
infinite lifetime (5.1 1). In this case the notion of iso-electronium 
essentially represents a tendency of pairs of valence electrons to 
correlate in singlet states at short distances. The use of iso- • 
chemistry, rather than conventional chemistry, is necessary 
because even this weaker form of the iso-electronium, as well as 
all screenings of the Coulomb potential at large, are non-unitary 
images of the conventional Coulomb settings as established in 
Section 5. As such, all these models require a representation on 
iso-spaees over iso-fields for their invariant formulation. In the 
text we present the projection of such an invariant formulation 
on conventional spaces over conventional fields for simplicity. 



7. Iso-chemical model of the hydrogen molecule with 
unstable iso-electronium 

In this section we study the solution of the restricted 
iso-chemical model of the hydrogen molecule, eqn (6.3), 
via conventional variational methods [1-4]- For this pur- 
pose we note that the solution of the full model with the 
Hulten potential exp( — rb)/[\ — exp( — rb)], where 
r, = b~ 1 is the horizon, implies rather considerable tech- 
nical difficulties. Therefore, we shall study model (6.3) 
under an approximation of the Hulten potential given by 
one Gaussian of the type [1 — A exp(— br)]/r with A a 
constant identified below. 

Recall from Section 5 that the stable character of the 
iso-electronium. eqn (5.1 1). is crucially dependent on the 
use of the attractive Hulten potential which 'absorbs' 
repulsive Coulomb forces at short distances resulting in a 
strong attraction. Therefore, the weakening of the Hulten 
potential into the above Gaussian form has the direct 
consequence of turning the iso-electronium into an 
unstable state. 

In this and in the following sections, we shall, therefore, 
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study a new iso-chemical model of the hydrogen molecule 
which is somewhat intermediary between the con- 
ventional chemical model and the iso-chemical model 
with a fully stable iso-electronium. 

It should be indicated that the assumption of an 
unstable iso-electronium implies the capability for more 
hydrogen atoms to bond to the H-molecules, as well as 
other features (such as the capability for the hydrogen 
and water molecules to acquire a paramagnetic character 
under external magnetic fields) for values depending on 
the lifetime of the iso-electronium. These features will be 
studied in detail in subsequent papers. 

The main objective of this section is to show the 
achievement of the exact representation of molecular 
characteristics even for the case of one Gaussian approxi- 
mation of the Hulten potential. Since a sufficient number 
of Gaussians can reproduce any potential, the achieve- 
ment of the exact representation with one Gaussian 
assures its preservation for additional Gaussians, thus 
also for the Hulten potential. 

The question whether the iso-electronium is stable or 
unstable evidently depends on the amount of instability 
and its confrontation with experimental data, e.g., on 
magnetic dipole moments. As such, the issue will be 
addressed theoretically and experimentally in a future 
paper. 

We begin with the approximation of the Hulten poten- 
tial with the Coulomb one, and then add a Gaussian 
'screening'. The reader should be aware that the latter 
brings us squarely within the applicability of iso-chem- 
istry because the map from the Coulomb potential to 
Gaussian screened potentials is generally characterized 
by a non-unitary transform, thus exiting the class of 
equivalence of quantum mechanics and chemistry. In 
fact, model (6.3) was constructed precisely via a non- 
unitary transform of the conventional quantum model. 

Under the above assumption, our first step is, there- 
fore, the study of model (6.3) in an exemplified Coulomb 
form characterized by the following equation, hereon 
expressed in atomic units (au) (see Ref. [21] for details): 

H x ¥ = ( - V 2 /4 - 2/r a - 2\r h + 1 \R) x (7. 1 ) 

where the differences from the corresponding equation 
for the Ht ion [1-4] are the replacement of the reduced 
mass (fi = ) 2 with 4, and the increase of the electric charge 
from (e = ) 1 to 2. 

The standard method for solving the above equation 
is the following [1-4, 32, 33]. The variational calculation 
is set up in matrix algebra form in a non-orthogonal basis 
set which has been normalized to 1. The metric of this 
non-orthogonal system of equations {S) is used to set up 
the orthogonal eigenvalue problem and the eigenvalues 
are sorted to find the lowest value. H is a Hermitean 
(square) matrix as well as Cand S; £is a diagonal matrix 
with the energy eigenvalues 
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HC = ESC; define C=S l 2 C\HS^ ' 2 C = ES ii2 C\ 
(S- l2 HS +i,2 C) = HC = E{S- in SS~ il2 )C = EC\ 

(7.2) 

where the last equation is obtained by multiplying the 
first equation from the left by 5" 1 2 and using the unitary 
property that S~ ' 2 = 5 + 1 2 to form an orthogonal eig- 
envalue problem. Finally we solve C by diagonalizing 
H and obtain C = S~ ,2 C. 

The detailed calculations have been done for both the 
restricted three-body model of the hydrogen molecule, as 
well as the full four-body versions. They are rather leng- 
thy and cannot be presented here for brevity. We shall, 
therefore, only outline the main results below and refer to 
[22] for a comprehensive presentation, including related 
references. 



8. Summary of results 

In order to demonstrate the advantage of the iso- 
chemical model using a Gaussian-screened-Coulomb 
attraction between electrons, a standard Boys-Reeves 
[22] calculation was carried out. This included all single- 
and double-excitations (CISD) from the ground state 
Hartree-Fock-Roothaan SCF orbitals for a 99x99 
'codetor' [23-28] interaction. Only the Is orbitals were 
optimized with a scaling of 1 . 1 9 1 for the least-energy 6G- 
ls orbitals, but the basis also included lG-2s, 2G-2p, 
lG-3s, lG-3p, 3G-3d, TdlG-4sp (tetrahedral array of 
four Gaussian spheres), and 4G-4f orbitals scaled to 
hydrogenic values as previously optimized. 

The additional basis functions provide opportunity to 
excite electrons to higher orbitals as is the standard tech- 
nique in configuration interaction (contrary to the main 
hypothesis of this work which is that there is an attractive 
'hadronic' force for electrons inside the r r critical radius). 
The results of the above calculations are summarized in 
Table 1 below. 

The Boys-Reeves C.I. achieved an energy of 

— 1.14241305 hartrees based on an SCF energy of 
-1.12822497 hartrees. This was followed by one 
additional iteration of 'natural orbitals' (C1NO) in which 
the first-order density matrix is diagonalized to improve 
the electron pairing to first-order [29]. The fact that this 
procedure lowered the energy only slightly to 

— 1.14241312 hartrees (-7.0x10' hartrees) indicates the 
99-configuration representation is close to the lower 
energy bound using this basis set while the iso-chemistry 
calculation produced the same exact energy with a com- 
paratively much smaller basis set. 

Since the Santilli-Animalu-Shillady LOBE (SASLOBE) 
has only an n 1 routine for the necessary integral trans- 
formation instead of the most efficient n s algorithm (w is 
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Table 1 

Summary of results for the hydrogen molecule 



Species H ; H ; H 2 

basis screening 



1 i- 

Is 




6.103 


1.191 


2s 


0.50 


24.35 


0.50 


2p 


0.50 


24.35 


2.36 


3s 


0.34 


16.23 




3p 


0.34 


16.23 




3d 


0.34 


-16.2' 




4sp 


0.25 


12.18 




4f 


0.25 


12.18 




Variational energy (au) 


* 


-7.61509174 




SCF energy (au) 


-1.12822497 


* 


- 1.13291228 


CI energy (au) 


-1.14231305 


• 




CINO energy (au) 


-1.14241312 






SAS energy (au) 


* 


* 


-1.174444 


Exact energy (au) [30] 


-1.174474 


* 


-1.174474 


Bond length (bohr) 


1.4011 


0.2592 


1.4011 


Electronium radius (bohr) 


* 




0.01124995 



1 The negative 3d scaling indicates five equivalent three-sphere scaled to 16.20 rather than 'canonical' 3d shapes. 

2 Three-body Hamiltonian (5.1). 



the number of basis functions), the SASLOBE C.I. runs 
are somewhat slow and required about 20 h on a 300 
M FLOPS Silicon Graphics computer. 

With more efficient routines this time can be reduced 
to about three hours. However, the screened -Coulomb 
attraction method used a smaller basis and achieved 
lower energies in a few seconds. It is also estimated that 
careful spacing of fewer quadrature points in the new 
integral routine can certainly reduce the SASLOBE run 
times by a factor of 2 at least. 

Therefore, it is clear that iso-chemistry calculations are 
(conservatively) at least 1000 times faster than a C.I. 
calculation. Another estimate is that since the new inte- 
gral corrections require a little more time than the usual 
Coulomb integrals (but do not take any additional stor- 
age space), the computer run-times for an iso-chemistry 
calculation should only be about three times the run- 
times for the corresponding Hartree-Fock-Roothaan 
calculation in any given basis set. 

The extension of the iso-chemical model of the H 2 
molecule to other molecules is conceptually straight- 
forward. In particular, the notion of iso-electronium 
essentially restricts all possible bonds to the established 
ones, as it is the case for the water molecule. 

In order to generalize the underlying quantitative treat- 
ment to molecules containing H-to-F, the pole-spike was 
re-optimized to obtain 100% of the correlation energy 
below the SCF energy in the given basis set since the SCF 
energy here was not quite at the Hartree-Fock limit (see 
Table 1 for details). 

The energy obtained here results from the calibration 



of the pole-spike to the experimental value of HF and is 
below the Quantum Monte Carlo energy of Luchow and 
Anderson [33] which requires hours on a much larger 
computer, as compared to less than 10 min for this work. 
In fact, the run times for HF were abut 8 cpu min on a 
100 MFLOP Silicon Graphics fc 02* workstation. 

The principal value of the pole (9.5 d) was calibrated 
for 100% energy of HF, H 2 0 has two tight sigma bonds 
and two diffuse lone-pairs so a single compromise value 
is a good test of the method. In HF the F" is nearly 
spherical so an average r r value does a better job of 
describing the 'correlation hole of transient iso-elec- 
tronium. 

The computed dipole moments are in excellent agree- 
ment with the experimental values. The use of the same 
pole value for H.O and HF degrades the H : energy 
slightly. The results of our studies for H-to-F based mol- 
ecules are summarized in Table 2 below. 

A comparison of the above data (particularly those on 
computer times) with corresponding data obtained via 
conventional approaches is instructive. 

9. Concluding remarks 

The fundamental notion of the new model of molecular 
bonds introduced in this paper is the bonding at short 
distances of pairs of valence electrons from two different 
atoms into a singlet quasi-particle state we have called 
iso-electronium, which travels as an individual particle 
on an oo-shaped orbit around the two respective nuclei. 
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Table 2 

Iso-electronium results for selected molecules 





H 2 


H 2 0 


HF 


SCF-energy (DH) 


-1.132800' 


-76051524 


-100.057186 


Hartree-Fock 4 






-1 00.071 85 4 


Iso-energy 


— 1.17444I 3 


-76J98229 3 


-100-459500' 


Horizon /UA) 


0.00671 


0.00038 


0.00030 


QMC energy 4 5 


-1.17447 


-76.430020 5 


-100.44296 4 


Exact non-rel. 


-1.174474 6 




- 100-4595 4 


% correlation 


99.9 2 


91 .6 2 


103.8 


SCF-dipole(D) 


0.0 


1-996828 


1.946698 


Iso-dipole (D) 


0.0 


1.847437 


1.841378 


Exp. dipole 


0.0 


1.85 7 


1.82 7 


Times (min : s) 


0:15.49 


10:08.31 


6:28.48 



(HD + )Dunmng-Huzinaga (10S/6P), [6,2, 1, 1, 1/4, 1. I] + H : P, + 3D1. 
' LEAO-6G1S + optimized GLO-2S and GLO-2P. 

2 Relative to the basis set used here, not quite HF-limit. 

3 Iso-energy calibrated to give exact energy for HF. - 

4 Hartree-Fock and QMC energies from Luchow and Anderson. 
5 QMC energies from Hammond et al. [22]. 

6 First 7 sig. fig. from Kolis and Wolniewicz [22]. 

7 Data from Chemical Rubber Handbook, 61st ed., p. E60. 
Horizon radius in Angstroms. 

Run limes on an 02 Silicon Graphics workstation (100 MFLOPS max.). 



The iso-electronium and related methodology are then 
assumed to characterize a covering of contemporary 
chemistry we have called iso-chemistry. 

The attractive short-range interactions needed to over- 
come the repulsive Coulomb force in the iso-electronium 
structure originate from non-linear, non-local and non- 
Hamiltonian effects in deep wave-overlappings; they are 
described by the recently achieved covering of quantum 
mechanics known as hadronic mechanics [1 1J; and their 
invariant formulation is permitted by the recently 
achieved broadening of conventional mathematics called 
iso-mathematics [10]. 

Specific experimental studies are needed to confirm the 
existence of the iso-electronium, by keeping in mind that 
the state may not be stable outside a molecule in which 
the nuclear attraction terms bring the electron density to 
some critical threshold for binding. 

Non-relativistic studies yield a radius of the iso-elec- 
tronium of 0.69 x 10" 10 cm. This 'horizon' is particularly 
important for iso-chemical applications and devel- 
opments because outside the horizon the electrons repel 
one another, while inside the horizon there is a 'hadronic 
attraction'. The value of the hadronic horizon obtained 
in this paper for H 2 via a Gaussian-lobe basis set (0.00671 
A) confirms the value independently obtained via non- 
relativistic studies. Nevertheless, the value of the had- 
ronic horizon must also be resolved via experiments- 
Similar considerations apply also for the depth of the 



hadronic well of the iso-electronium which has been esti- 
mated at about 1 MeV. 

The same non-relativistic studies also predict that the 
iso-electronium is stable within a molecule, although par- 
tially stable configurations also yield acceptable results. 
The question of the stability vs instability of the iso- 
electronium must, therefore, also be left to experimental 
resolutions. 

The foundations of the iso-electronium can be seen in 
a paper by Santilli [29] of 1978 on the structure of the 
7f-meson which contains the first identification of the 
attractive character of non-linear, non-local and non- 
Hamiltonian interactions due to deep wave-overlappings 
in singlet coupling (and their repulsive character in triplet 
coupling). The iso-electronium also sees its foundations 
in subsequent studies by Animalu [30] of 1994 and Ani- 
malu and Santilli [31] of 1995 for the Cooper pair in 
superconductivity, as well as in other examples of electron 
bondings existing in nature, such as ball lightning. 

The iso-electronium also results in having deep con- 
nections with a variety of studies in chemistry conducted 
throughout this century [23-28, 32, 33], and actually pro- 
vides the physical-chemical foundations with a more 
adequate mathematical formulation for most of them. 

Therefore, the iso-chemical model of molecular bonds 
results in being consistent with a number of important 
occurrences in particle physics, superconductivity, chem- 
istry and other fields, such as: 
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( 1 ) The ro-;iectroniuni introduces a new attractive force 
among the neutral atoms of a molecular structure 
which is absent in quantum chemistry and certainly 
improves our understanding of the strength of molec- 
ular bonds and its stability. 

(2) The iso-electronium permits an immediate interpret- 
ation of the reasons why the H : and H : 0 molecules 
only admit two H-atoms. 

(.V) The iso-electronium permits tiie achievement of a 
representation of the missing 2% in the binding 
energy, thus allowing meaningful thermodynamical 
calculations. 

(4) The iso-electronium provides an explanation of the 
long known, yet little understood Pauli exclusion 
principle, according to which electrons correlate 
themselves in singlet when on the same orbit without 
any exchange of energy, thus via interactions essen- 
tially outside the representational capabilities of 
quantum mechanics and chemistry. 

(5) The tso-electronium is consistent with the known - 
existence of superconducting electron-pairs which 
bor.d themselves so strongly to tunnel together 
through a potential barrier. 

(6) The iso-electronium provides a quantitative model 
for the explanation of "electron correlation*. Instead 
of a complicated 'dance of electrons" described by 
posithe energy excitations, the iso-chemistry expla- 
nation is that electrons are energetically just outside 
the horizon of a deep attractive potential well due to 
their wavefunctions overlapping beyond the critical 
threshold of the hadronic horizon. 

(7) The iso-electronium is consistent with the 'Coulomb 
hole' documented by Boyd and Yee as found from 
subtracting accurate explicitly -correlated wav- 
efunctions from self-consistent-field wavefunctions. 
In our studies the 'hole' is re-interpreted as a 'had- 
ronic well 1 . 

(8) The iso-electronium is also in agreement with the fc bi- 
polaron* calculated for anion vacancies in KG by 
Fois. Sellonu Parinello and Car and bi-polaron spec- 
tra reported by Xia and Bloomfield. 

(9) The iso-electronium permits an increase of the speed 
in computer calculations at least 1 000-fold. 

As we hope to illustrate m future works, a promising 
feature of the proposed iso-chemical model of the hydro- 
gen molecule is not only the capability to permit accurate 
representations of experiment:-'.! data in shorter computer 
times, but also the capability to predict and quantitatively 
treat new features and reactions. 
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Abstract 



Despite outstanding advances throughout this century? we still lack final knowledge on the structure of the water 
molecule because of a number of insufficiencies in available models identified in the text. In a preceding paper we 
submitted a covering of quantum chemistry called hadronic chemistry in its isochemical branch, and introduce a 
new model of the hydrogen molecule characterized by a bond at short distances of the two valence electrons into a 
singlet quasi-particle state called isoelectronium. In this paper we introduce, apparently for the first time, the 
isochemical model of the water molecule which is characterized by two isoelectronia, one per each H-0 dimer and 
show that the model: (1) introduces a new strong force responsible for the molecular bond; (2) explains the reason 
why the water molecule has only two hydrogen atoms; (3) provides an accurate representation of the binding energy 
from the first axiomatic principles without ad hoc modifications; (4) permits a correct representation of electric and 
magnetic moments; and (5) reduces by at least a factor of 1000 computer usages for calculations due to a much 
faster convergent series. Experimental verifications, applications and novel predictions of the model are studied in 
subsequent papers. © 1999 International Association for Hydrogen Energy. Published by Elsevier Science Ltd. All 
rights reserved. 



1. Introduction 

Water is an extremely important compound on 
Planet Earth in a biological as well as geophysical 
sense. For this reason, it has been subjected to compre- 
hensive studies conducted since the beginning of quan- 
titative science with outstanding results (see, e.g., Refs. 
[1,4,14,38,17,21]). 

Despite that, a number of fundamental issues on the 
structure of the water molecule remain still open at the 
time of writing. Among them we note: the inability to 
explain why the hydrogen and water molecules only 
have two H -atoms (because the admitted forces, such 
as the exchange and van der Waals forces, are of 
nuclear type thus, being conceived and applicable for 
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an arbitrary number of constituents); the binding 
energy still misses 1 to 2% under the strict implemen- 
tation of the quantum axioms; accurate representations 
are based on structural changes of the Coulomb law, 
such as via Gaussian screenings whose axiomatic impli- 
cations require a specific study: larger percentages are 
missing in the representation of electric and magnetic 
moments (which sometimes even have the wrong sign); 
calculations continue to require large computer times 
because of poorly convergent series; we still miss an 
explanation of the electron 'correlation energy'; and 
other still open basic issues. 

When passing from the structure of one water mol- 
ecule to more general profiles, the number of open 
basic issues increases. For instance, it is generally 
admitted that quantum chemistry has been unable to 
provide a systematic theory of the liquid state in gen- 
eral, let alone that of the water-m particular. 

Also, chemical reactions"^ general are irreversible, 
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while the axiomatic structure of quantum chemistry is 
strictly reversible because the theory is strictly 
Hamiltonian and all known potential forces are revers- 
ible. This results in a seemingly irreconcilable incom- 
patibility of the very axiomatic structure of quantum 
chemistry with chemical reactions in general, and those 
involving water in particular. In fact, there is no 
known rigorous way in which the reversible axiomatic 
structure of quantum chemistry can be turned into the 
irreversible form needed to represent irreversible events 
of the chemical reality. 

When passing to water as a constituent of biological 
entities, the limitations of quantum chemistry reach 
their climax. In fact, biological structures such as a cell 
are not only irreversible (because they grow, age and 
die), but have such a complex structure to require 
multi-valued theories (also known in mathematics as 
hyperstructures). The expectation that quantum chem- 
istry, with its reversible and single-valued structure, 
can effectively represent biological structures and their 
evolution is beyond the boundaries of science. 

In view of the above limitations, in the preceding 
paper Ref. [2,24-29] we have constructed a covering of 
quantum mechanics under the name of hadronic mech- 
anics. In essence, the two mechanics coincide every- 
where, except the region inside a sphere of radius of 
the order of 1 fm (= 10~ 13 cm) called hadronic horizon 
in which the broad theory holds. Hadronic mechanics 
then results to be a form of 'completion' of quantum 
mechanics much along the lines of the historical 
Einstein-Podolsky-Rosen argument, although 
achieved via the addition of contact, nonHamiltonian, 
nonlinear, nonlocal and nonpotential forces due to 
deep overlapping of the wavepackets of particles. 

On more technical grounds, hadronic mechanics is 
based on new mathematics, called iso- geno- and 
hyper-mathematics for the characterization of general- 
ized reversible, irreversible and multivalued systems, 
respectively. 

On pragmatic grounds, hadronic mechanics can be 
easily constructed via a step-by-step transform UU^ of 
all aspects of quantum mechanics which is equal to the 
unit (unitary) everywhere, except inside the hadronic 
horizon (where it is nonunitary). NonHamiltonian 
effects due to deep wave overlappings can therefore, be 
easily added to a conventional model (when appropri- 
ate) via nonunitary transforms of the latter. 

Following the construction of hadronic mechanics, 
the authors have constructed a corresponding covering 
of quantum chemistry under the name of hadronic 
chemistry [3,30], which also comprises three different 
branches, isochemistry for the representation of revers- 
ible yet generalized structures; genochemistry for the 
representation of irreversible structures; and hyperch- 
emistry for the representation of irreversible multi- 
valued structures. 



Isochemistry was then applied for the construction 
of a new model of the hydrogen molecule and proved 
that it does resolve at least some of the insufficiencies 
of quantum chemistry, such as an explanation that the 
molecule has only two H-atoms; a representation of 
the binding energy exact to the seventh digit; computer 
calculations which converge at least 1000 times faster 
than those of quantum chemistry and other advances. 

In this paper we shall study the isochemical branch 
of hadronic chemistry; we shall submit, apparently for 
the first time, the isochemical model of the water mol- 
ecule; and we shall show that it does indeed permit the 
resolution of at least some of the insufficiencies indi- 
cated above. 

A technical understanding of this paper, requires at 
least some knowledge of Refs. [2,24-29,3,30] whose 
review is omitted here for brevity. 



2. Exactly solvable isochemical model of the water 
molecule with stable isoelectronium 

In the preceding paper Ref. [3] we have introduced 
the main hypothesis of the isochemical molecular 
model according to which two electrons from two 
different atoms bond themselves into a singlet quasista- 
ble and quasiparticle state called isoelectronium which 
describes an oo-shaped orbit around the nuclei, as it is 
the case for planets in certain binary starts. 

The main characteristics of the isoelectronium in the 
first nonrelativistic approximation were calculated and 
resulted in being: charge -2e, spin 0, magnetic dipole 
moment 0, mass 1.022 MeV and radius 

r c = b~ l *{h 2 lmV) Xil = (h/ma> 0 ? /2 

= 6.8432329 x 10" n cm = 0.015424288 bohrs 

= 0.006843 A. (1) 

In the above nonrelativistic approximation, the mean 
life resulted in being infinite (full stability, see Ref. [3], 
Section 5 in general and Eq. (5.11) in particular), with 
the understanding that relativistic corrections are 
expected to render it finite (partial stability). AH con- 
ventional forces of current use in chemistry (the 
exchange, van der Waals and other forces) then hold 
between the isoelectronium and the remaining elec- 
trons. 

The model was applied to the isochemical model of 
the hydrogen molecule and implied: (1) the introduc- 
tion of a new bonding force given by the Hulten force 
besides those in current use, which is so strong to over- 
come the repulsive Coulomb force (see Ref. [3], 
Section 5 in general and Eq. (5.6) in particular), thus 
being capable of representing the strong molecular 
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Fig. 1. A schematic view of the proposed isochemical model 
of the water molecule which is depicted in terms of orbits of 
the valence electrons (rather than in terms of density distri- 
butions) when at absolute zero degrees temperature and in the 
absence of rotations, vibrations, translations and other 
motions. The fundamental assumption is that the two valence 
electrons, one per each pair of atoms, correlate themselves 
into a bonded singlet state at a short distance we have called 
isoelectronium, which is assumed to be stable. In this case the 
only orbit yielding a stable water molecule is that in which 
the isoelectronium describes an oo-shaped orbit around the 
respective two nuclei of the H- and O-atoms. The isoelectro- 
nium is then responsible for the attractive force between the 
two atoms. The binding energy is instead characterized by the 
oo-shaped orbit of the isoelectronium around the two nuclei, 
conceptually represented in this figure via a standing wave for 
a particle of spin 0 and charge -2e. In this way the water 
molecule result in being characterized by two isoelectronia, 
one per each dimer H-O. The water molecule is then reduci- 
ble to two intersecting H-O dimers with a common O-atom. 
Note that, already in the absence of molecular motions, the 
orbits of the two isoelectronia are perpendicular to the H-O- 
H plane, thus confirming a characteristic of the water mol- 
ecule indicated in Ref.[l]. Needless to say, when rotations are 
added, they imply the rotation of the orbits of the isoelectro- 
nia around their individual symmetry axis, resulting in two 
spherical distributions of the electrons, thus recovering con- 
ventional spherical charge distributions. Conventional 
exchange, van der Waals and other forces then imply the 
reconstruction of the conventional charge distributions [}]. 
The model then permits a representation of: (1) the strong 
value of the molecular bond; (2) the reason why the H 2 0) 
molecule has only two hydrogen atoms; (3) the exact rep- 
resentation of the binding energy, electric and magnetic 
moments; and other advances studied in the text. The model 
of the H-O dimer studied in the text is then extendable to 
other H -based dimers, such as H-C 



bonds as occurring in reality; (2) the immediate expla- 
nation of the reason why the molecule has only two 
hydrogen atoms, evidently because, once the two elec- 
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trons of the H 2 molecule are bonded into the isoelec- 
tronium, there is no possibility of bonding a third one; 
(3) a representation of the binding energy of the 
hydrogen molecule accurate to the seventh digit; (4) a 
decrease of computer usage for calculation at least 
1000 fold; and (5) other advances studied in sub- 
sequent papers, such as the impossibility for the mol- 
ecule to acquire a net magnetic polarity because the 
two o-branches of the isoelectronium orbit imply 
opposite rotational motion, with consequently opposite 
magnetic polarities in agreement with reality. 

In this paper we introduce, apparently for the first 
time, the isochemical model of the water molecule 
H 2 0 = H-O-H considered at °C and in the absence of 
any rotational, oscillation or other motion. The main 
hypothesis is that each electron from the two H-atoms 
couples in singlet with one 2p electron from the O- 
atom, resulting in two isoelectronia, one per each H-O 
dimer as in Fig. 1 . 

According to this model, the bonding force between 
the two H- and O-atorns is given by the Hulten force 
between the two valence electrons in the isoelectronia, 
plus contributions from the conventional exchange, 
van der Waals and other forces. The binding energy is 
instead characterized by the two oo-shaped orbits of 
the isoelectronia around the H-O-H nuclei. 

Again, in the first nonrelativistic approximation we 
shall assume in this section (see next section for a 
different assumption) that the two isoelectronia of the 
water molecule are permanently stable and have 
charges -2e, spin 0, magnetic dipole moment 0, and 
the same radius as that in the H 2 molecule, Eq. (1). On 
the contrary, the radii of the two oo-branches of the 
H-O dimers are expected to be different between them- 
selves, as well as different than the corresponding 
value in the H 2 molecule. 

In this section we shall introduce, apparently for the 
first time, a hadronic/isoschroedinger equation for the 
water molecule which is evidently approximate, yet 
exactly solvable. We shall then show that the model is 
extendable to all other dimers comprising one hydro- 
gen atom, such as H-C. 

For this purpose, we approximate the H-O-H mol- 
ecule as being composed of two intersecting identical 
dimers H-O with evidently only one oxygen atom. 
This requires a first correction due to the lack of inde- 
pendence of said dimers reviewed in Section 3. 
Moreover, in each H-O dimer we shall assume that 
the oxygen appears to the isoelectronium as having 
only one net positive charge + e located in the nucleus. 
This evidently requires a second correction which 
essentially represents the screening of the various elec- 
trons of the oxygen. Additional corrections are also in 
order along conventional lines [1]. 

A study of these corrections has indicated that they 
can all be represented via one single Gaussian screen- 
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ing of the Coulomb law of the type 
+e +e(l±e-°^) 



(2) 



where p is a positive parameter to be determined from 
experimental data, the sign - applies for the screened 
O-nucleus as seen from an electron (because of the 
repulsion caused by the electron clouds of the oxygen), 
while the sign + applies for the screen O-nucleus as 
seen from the H-nucleus (because of the attraction this 
time caused by said electron clouds). 

The resulting model is structurally equivalent to the 
isochemical model of the hydrogen molecule of Ref. 
[3], except for the modifications indicated about, and 
can be outlined as follows. 

By denoting with the sub-indices 1 and a the hydro- 
gen, 2 and b the oxygen, prior to the indicated screen- 
ing and in the absence of all hadronic effects, the 
conventional Schroedinger equation of the H-0 dimer 
with the oxygen assumed to have only one elementary 
charge + e in the nucleus is given by 
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As it was the case for the H 2 -molecule, our task is that 
of subjecting the above model to a transform which is 
nonunitary only at the short mutual distances 
r c = b~ x = r n of the two valence electrons (here 
assumed to be inside the hadronic horizon), and 
becomes unitary at bigger distances I r g 10" 13 cm / J, 
I r ^>10~ 13 cm = /. This guarantees that our isochemical 
model coincides with the conventional model every- 
where except for small contributions at small distances. 

We assume that the state and related Hilbert space 
of systems (3) can be factorized in the familiar form 
(in which each term is duly symmetrized or antisymme- 
trized) 

I =1 ^i2>x I *i.)x I *ib>+ I I ^2b>* I *r>> 

jr Tox = Jf'n x 3f u x Jf ]b x ^2a x Jf 2 b x jf R . (4) 

The nonunitary transform we are looking for shall act 
only on the r n variable characterizing the isoelectro- 
nium, while leaving all other variables unchanged. The 
simplest possible solution is given by 

l/(^ I2 )UVl2) = I = «?* (,te)/^<n,),^i,n * ,<hl),1 "* ( ' d,,, . (5) 

where the ^'s represent conventional wavefunction and 
the represent isowavefunctions. 



We now construct the isochemical model by trans- 
forming short-range terms (isochemistry) and adding 
untransformed long-range ones (chemistry), thus 
resulting in the radial equation 
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By using the derivation of Ref. [3] and by recalling 
that the Hulten potential behaves at small distances 
like the Coulomb one, Eq. (6) becomes 
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The above model can be subjected to an important 
simplification. In fact, in the first approximation herein 
considered, the H-O dimer (7) can be reduced to a 
restricted three body problem similar to that possible 
for the conventional H 2 + ion [1], but not for the con- 
ventional H 2 molecule, according to the equation 
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The indicated corrections due to the screening of the 
various electrons in the oxygen and other corrections 
are needed in the 'sensing' of the O-nucleus by the iso- 
electronium as well as by the H-nucleus, yielding in 
this way our final model 
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(9) 



where a is a positive parameter and E is half of the 
binding energy of the water molecule. 

Under the latter approximation, the model permits, 
for the first time to our knowledge, the achievement of 
an exact solution for the structure of the water mol- 
ecule which however, exists only for the case of the 
fully stable isoelectronium. In fact, for the unstable 
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isoelectronium the model becomes a four-body struc- 
ture which, as such, admits no exact solution. 

Besides being exactly solvable, model (9) exhibits a 
new explicitly attractive force among the neutral atoms 
of the H-O dimer which is absent in conventional 
quantum chemistry; the equation also explains the 
reasons why the water molecule admits only two H- 
atoms; as we shall see in the remaining sections, the 
model also permits the exact representation of the 
binding energy, electric and magnetic moments; the 
model also yields much faster convergence of a series 
with much reduced computer times, and resolves other 
insufficiencies of conventional models. 

Finally, the model is evidently extendable with 
simple adjustment to an exact solution of other dimers 
involving the hydrogen, such as H-C. In addition, it 
permits the identification of electric and magnetic 
polarizations which are not predictable with quantum 
chemistry. 



3. Gaussian approximation of the isochemical model of 
the water molecule with unstable isoelectronium 

The solution of the exactly solvable Eq. (9) will be 
studied elsewhere. In this section we shall study a 
Gaussian approximation of the isochemical model of 
the water molecule in which the Hulten potential is 
essentially approximated to a certain Gaussian form. 

It should be indicated from the outset that such an 
approximation implies an evident weakening of the 
Hulten attraction among the two isoelectrons of the 
isoelectronium which, in turn, implies the instability of 
the isoelectronium itself, thus reaching a model which 
is somewhat intermediate between the full isochemical 
model and the conventional quantum mechanical 
model of the water. 

Despite this approximate character, the results of 
this section are significant because they show the capa- 
bility by isochemistry to achieve an essentially exact 
representation of the binding energy, electric and mag- 
netic moments and other characteristics of the water 
molecule. 

But any potential can be exactly reproduced by a 
sufficient number of Gaussians, including the Hulten 
potential. Therefore, the exact representations of this 
section reached for only one Gaussian can be evidently 
expected to persist for the full isochemical model. 

Finally, Gaussian and other screenings of the con- 
ventional Coulomb law are now known to be outside 
the equivalence class of quantum chemistry because 
they require a nonunitary transform [3]. The Gaussian 
approximation of this section is, therefore, fully within 
the basic, nonunitary, axiomatic structure of hadronic 
chemistry, while similar screenings have questionable 



axiomatic foundations when believed to be part of 
quantum chemistry. 

The results of this second study can be outlined as 
follows. Since HOOH will be slightly more 'allowed' 
under collisions of HO" with neutral H 2 0 and the in- 
ternal repulsion within the anion could favor the 
release of (2e)" 2 to form OH + . Collisions of OH + 
with OH" will then further enhance the concentration 
of HOOH and transport of (-2e) will contribute to 
the current. 

The question here is whether under extreme cases of 
forced conduction a singlet-pair of electrons (isoelec- 
tronium, can be 'triggered' (in the nomenclature of 
Santilli [6] and Animalu and Santilli [31]) within a 
water molecule to form and release a -2e charged iso- 
electronium particle which will provide an additional 
conduction mechanism analogous to Cooper-pairs of 
electrons in superconducting solids. Since the energy 
depth of the V 0 parameter in the isoelectronium 
Hulten potential of the original 1978 derivation by 
Santilli [6] and Animalu and Santilli [31] is not known, 
nor how closely the Gaussian representation fits the 
Hulten form, we can only match the radius of the two 
potentials and calculate the energy differences caused 
by the 'sticky-electrons' model in which a transient 
form of isoelectronium can occur (the Gaussian poten- 
tial may well not be deep enough to ensure a perma- 
nent bound state for isoelectronium). 

The 'sticky-electron' model is a parametric model 
which includes the magnetic dipole attraction 
between singlet-paired electrons as well as the nonlo- 
cal merging of the wave-packets of each electron at 
short distance. The radius of the Gaussian screening 
is then determined empirically by fitting the calcu- 
lated energy as nearly as possible to the most accu- 
rate energy values available. As used here it should 
be emphasized that the off-axis positions of the 
Gaussian-lobe basis sets [7,8,32-36,9,43,37] ensure 
that angular correlation is included as well as radial 
dependence and can include the magnetic dipole 
attraction of opposite electron spins as well as mer- 
ger of wavepackets. One radial screening parameter 
used with off-axis basis sets parametrically covers all 
forms of short range attraction which may include 
angular dependence. Thus, the present model can 
give us an approximate energy difference required to 
release an electron-pair from OH". 



OH"— >OH + - 2e. 



(10) 



It will be seen below that the energy difference between 
OH" and OH + as calculated, allowing a transient 
form of isoelectronium, is well within the voltages 
accessible using capacitive discharge through water. 
Such a mechanism which would allow (-2e) particles 
to flow through water would not be superconductivity 
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since the freely moving molecules and ions are not 
constrained to lattice positions as in solids so that 
resistive I 2 R heating will still occur. This is mainly due 
to the fact that conduction in liquids occurs by mobi- 
lity of both anions and cations along with size differ- 
ences, polarizability differences and special mechanisms 
such as the hydrogen -bonding 'domino' effect for H + 
transport. In solid-state conduction only the electrons 
move (ignoring in-place phonon oscillations because 
the atoms do not travel from one electrode to the 
other; the apparent motion of positively charged 
'holes' is also due to the motion of electrons while the 
atoms merely oscillate about the mean lattice pos- 
itions), but in solutions there is two-way traffic with 
positive and negative ions traveling in the opposite 
directions and with differing velocities; this will cer- 
tainly lead to resistive heating even up to the vaporiza- 
tion of the water as well as a high probability of ion - 
collisions. It should be noted that in recent work 
Zhitenev et al. [10] have measured migration of paired- 
electrons to quantum dot wells in Gaas and Boyd and 
Yee [11] have observed 'bipolaron* electron pairs in 
alkali halide lattice vacancies. Calculations leading to 
unexpected bipolarons in crystal lattice vacancies have 
also been observed by using the method of Fois et al. 
[12]. These findings in solids lend support to the con- 
cept of an electron pair as an individual particle, called 
isoelectronium here. 

The calculations given here do not prove the pre- 
sence of isoelectronium particles in high current aqu- 
eous electrical conduction, they only indicate the 
energy threshold necessary to form (-2e) within the 
conducting solution by double-ionization of OH~. It is 
not easy to envision an experiment that would be able 
to analyze components of a given current due to a 
multiple ion species in terms of the amount of current 
due to (— 2e) particles and none is proposed here. 
However, there may be a chemical test for such a 
mechanism. Once OH~ is double-ionized to form 
OH + , collisions with (— 2e) particles would regenerate 
OH~ ions just as collisions of H + with OH~ will 
reform H 2 0 and no new species will be evident. 
However, if OH~ collides with OH + a new chemical 
species (HO-OH) will be formed that may last long 
enough in the 'liquid plasma' to behave as a strong 
oxidizing agent. Thus, organic compounds with double 
bonds (alkenes) which have negligible conductance 
could be added to water undergoing a high current 
flow to cause hydroxylation of such compounds [13] 
(i.e., conversion of alkenes to epoxides which are then 
readily hydrolyzed in the presence of H + to diols). 
Enhanced concentrations of epoxides and diols would 
be indirect evidence of double ionization of OH~. 



OH* + OH~-» HO-OH 
CH 2 =CHR+ HO-OH -»CH=CHR + H,0 

O 

(11) 

If isoelectronium can be detected indirectly by a chemi- 
cal method this would in itself be an important infer- 
ence on the existence of a two-electron, spin -paired 
particle. More important, 'isochemical reactions* could 
be driven by high conduction 'liquid plasma' environ- 
ments where isoelectronium is at an enhanced concen- 
tration. 

Another case of interest is that of aqueous 'mixtures' 
of insoluble organic compounds forming a separate 
'oil layer' over water in an intense magnetic field of 
several Tesla. At normal thermal energy of room tem- 
perature kT (RT per mole) the main energy form 
would be a random Brownian motion. However, in the 
presence of a strong magnetic field HO* and HO~ 
would be constrained to favor circular motions in the 
magnetic field by the 'cyclotron effect*, but there is no 
obvious source of HO + . Since two ions of opposite 
sign charges would be favored to collide by both elec- 
trical attraction and by opposite path curvature in a 
magnetic field, there is an enhancement that when cre- 
ated as a normal result of the (H + , OH", H~, H 2 0 2 , 
HO J, H 2 0) equilibrium system studied by Pourbaix 
[5], any natural concentration of HOOH would be 
augmented by collision of H + with HOO". In addition 
there is some slight change that H + would collide with 
OH" with sufficient excess energy to produce OH + 
and H". Thus, the presence of an intense magnetic 
field 

H + + HOO--^HOOH, (12a) 
H + -f OH" — H" + OH+ (1 2b) 

OH + + OH" -►HOOH (12c) 

would cause positive and negative ions to collide more 
easily while traveling in the opposite rotational arcs in 
such a way as to enhance the concentration of HOOH 
which could then epoxidate alkenes and upon hydroly- 
sis would lead to diols. A direct measure of this effect 
would be to determine the enhanced solubility of 
alkenes in water. The alkenes are only slightly soluble 
in water ('oil and water do not mix*) but alkenes con- 
verted to diols will have a measurably greater solubility 
in water due to the attached -OH groups. Again, if 
such enhanced solubility of alkenes in water can be 
caused by intense magnetic fields, this would be indir- 
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ect evidence of the existence of an electron-pair par- 
ticle. 

Note below in the description of the calculations 
that the key to this work is that it is easy to calculate 
the energy of OH" when one subtracts a small amount 
from the two-electron repulsion terms in the usual 
HFR-SCF treatment due to the attraction of singlet- 
paired electrons at close range within 1.0 picometer. In 
the recent 'Handbook of Computational Quantum 
Chemistry* by Cook [14,38] (p. 438) it is noted that 
solutions to the HFR-SCF scheme may not always 
exist for anions! However, in the method used here 
convergence of the HFR-SCF process was normal for 
an SCF process, presumably because the so-called 
'self-energy' error of the Hartree-Fock method [15] 
(each electron repels all electrons including itself) is lar- 
gely cancelled by the new attractive terms used here. In 
effect this description of OH" is possible because of 
the easy convergence of the *correlated-SCF' process. 

The model adopted here is to use the usual Hartree- 
Fock-Roothaan self-consistent-field equations [16] 
(which also has some formal flaws such as the self-in- 
teraction terms [15]) and question the form of the 
Coulomb interaction of the electrons. Note that redu- 
cing the values of the Coulomb integrals will lower the 
energy by reducing the electron-electron repulsion 
while reducing the exchange terms will raise the 
energy, but the (1/2) factor reduces the effect of the 
exchange terms. Thus, a reduction of the value of the 
integrals will lower the energy. Note that Goodgame 
and Goddard [15] have already recommended reducing 
the atomic self-energy by subtracting 1.39 eV from 
Hartree-Fock exchange integrals in the cases of Cr 2 
and M0 2 . 

FC = ESQ 

F itj = H u + IXa, iWJ/k, /) (I3a) 
- 1/2(1, k/j, 1)1 

<W/*. /) = J j X,( I 1 )/( 1 Ai2)/Xj t (2)/ 1 (2)dr ! dr 2 , 

(13b) 

Pi tj = 2L n cn, ic nJ 

(13c) 

(sum n only over occupied orbitals). 

The 1995 paper on Cooper pairs by Animalu and 
Santilli [31] invokes the nonlocal hadronic attractive 
force evident in the 7t°-meson as applied to a pair of 
singlet-paired electrons which form a boson quasi-par- 
ticle. However, the 'collapsed* positronium rapidly 
decays since the particle-antiparticle annihilation takes 



place in less than a picosecond. In the electron-elec- 
tron case it is believed that there may be a stable 
quasi-particle. After using a nonlocal isotopic non- 
linear transformation, the hadronic attraction was 
transformed back to real-space and modeled with a 
Hulten potential. In this work considerable effort was 
made to evaluate the matrix elements for the Hulten 
potential without success. Examination of the original 
1978 paper on positronium collapse by Santilli [6] and 
Animalu and Santilli [31] revealed that the Hulten po- 
tential is not necessarily a unique representation of the 
hadronic force; at least a linear combination of similar 
potentials could be used to represent the Hulten if the 
matrix element of such other potentials could be evalu- 
ated. The depth of the screened Gaussian approxi- 
mation is determined here by requiring that the width 
at half height of the Gaussian is equal to the *b' value 
of the Hulten 'horizon', the radius at which the 
Coulomb repulsion is zero. Thus, the screened- 
Gaussian potential probably has a depth which is too 
shallow, although the V 0 depth parameter for the 
Hulten potential is not known at present. 

This work assumes that until the matrix elements of 
a two-electron interaction for singlet-pairs can be 
found for the Hulten potential, a Gaussian-screened- 
Coulomb potential can be used to describe the real- 
space form of the hadronic attraction and a parameter 
fitted to experimental energies the screening exponent 
probably includes other effects such as the magnetic 
dipole interaction of two electrons with opposite spin- 
magnetic-moments. This form has the important prop- 
erty that it can be merged with the general case of the 
four-center Coulomb or exchange integral derived by 
Shavitt [17] using the famous Gaussian transform tech- 
nique. 

The Gaussian transform two-electron integral for 
four Gaussian spheres has been used in a number of 
Gaussian-lobe basis SCF programs written by Mosier 
and Shiilady [18], Shiilady and Baldwin [39], Shiilady 
and Talley [19], Shiilady and Richardson [40] and 
others. It is important to note that the formula is com- 
pletely general in orientation of four Gaussian sphere 
lobe-orbitals as well as the distance between two elec- 
trons. As modified for description of correlation of 
two electrons, such a general formula can describe 
angular correlation as well as distance interaction. 
Thus, matrix elements of a screened-CouIomb inter- 
action were subtracted from the usual (1/r) Coulomb 
repulsion to model the real-space form of the hadronic 
attraction of two electrons. This work added the 
Gaussian screening as exp(— ar 2 )/r so that the special 
properties of Gaussians could be used, especially the 
properties that the product of two Gaussians form 
another Gaussian (times a rec-entering factor) and that 
polar coordinates readily separate into factorable x, y 
and z components. The goal is to evaluate the two- 
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Table 1 

Isoelectronium results for selected molecules 0 










OH + 


OH" 


H 2 0 


Hr 


SCF^energy (DH) 
Hartree-Fock c 
Iso-energy b 
Horizon RC (A) 
QMC energy^ 
Exact non-rel. 
Iso-dipole (D) 
Exp. dipole 


-74.860377 

-75.056678 
0.00038 

5.552581 


-75.396624 

-75.554299 
0.00038 

8.638473 


-76.058000 

-76.388340 

0.00038 
-76.430020 d 

1.847437 
1.84 [2] 


-100.060379 
-100.07185 c 
-100.448029 

0.00030 
-100.44296 c 
-100.4595 c 

1.8413778 

1.82 


a Notes: (HD + ) Dunning-Huzinaga (10 S/6P), [6,2,1,1,1/4,1,1] + H2S1 + H2P1 + 3DI. 

b Iso-energy calibrated to give maximum correlation for HF. 

c Hartree-Fock and QMC energies from Luchow and Anderson [22]. 

d QMC energies from Hammond et al. [21]. 



electron four-center matrix elements of the Gaussian- ~ 
screened-Coulomb potential shown below. 

Y(r) = (1-2 exp[-(a^)]/r. (14) 

Amazingly, the Gaussian-transform variable V can be 
simply shifted by the V value of the screened - 
Gaussian exponent and carried through the original 
derivation until the last step when integration over V 
is required. Alpha is usually a very high number, this 
work used 0.13441885 x .10 7 . At this point the usual 
Coulombic interaction resorts to a well-known auxili- 
ary function f 0 which has been studied by M osier and 
Shillady [18], Shillady and Baldwin [40] and others. 
Since both s 1/2 and (s+a) 1/2 occur in the denominator 
of the screened-Coulomb form, two poles occur in the 
integral. A change of variable absorbs the pole due to 
(5-f-<a)~ 1/2 and shifts the other pole due to (s)~ 1/2 to 
the lower limit of the integral A smooth spike is evi- 
dent at the lower limit of the numerical integration 
using a 70 point Simpson's rule integration (two ranges 
are used with 20 points more closely spaced near the 
pole and 50 points for the remaining range). 

This work was carried out using 64 bit double pre- 
cision arithmetic which provides 14 significant figures. 
A simple offset (<5) of 1.0 x 10~ 15 has provided useful 
results with this simple offset to avoid numerical over- 
flow. While this pole is formally a problem in needing 
a continuous function to integrate, numerical inte- 
gration seems to handle this well out to 14 significant 
figures; particularly since the routines used for the 
Coulombic integrals are known to be accurate only to 
12 significant figures. The area under the pole-spike is 
estimated as a narrow triangle upon a rectangle 
1.0 x 10" 15 wide with the height of the triangle set at 
3.43207 x 10 8 times the height of the point set 
1.0 x 10" 15 into the range of integration (the first 
Simpson point). The present code for this screened- 



Coulomb integral is presently slower than the corre- 
sponding foe [17] function used for the Coulombic 
integrals due to the 70 point Simpson integration, but 
the integrand is nearly flat after the spike at s = 0.0 so 
that the portion of the integrand can be evaluated 
more rapidly with fewer points. For the results pre- 
sented here, the simple offset of the lower limit by 
1.0 x 10" 15 is adequate for this report. Further details 
on the auxiliary integral can be found in a previous 
paper on the H 2 molecule [20,41,42]. Work in progress 
indicates it may be possible to express the new auxili- 
ary integral to an analytical expression involving the 
Erf(x) function, but until further checks are completed 
this work used the Simpson integration. Note the inte- 
gral is a result of a simplification of a 12-fold inte- 
gration over the volume elements of two electrons and 
has been reduced to a one-dimensional integration 
multiplied by appropriate factors. 



4. Results 

The geometry given for H 2 0 by Dunning [20,41,42] 
was used to carry out the usual HFR-SCF calculation 
after an additional 3 d orbital mimic [19,40] was opti- 
mized for the O-atom and (2 S,2p) orbitals were added 
for the H-atoms. The exponent for the 03 d orbitals 
was optimized to three significant figures and the 
(03 d,H2 s,H2p) exponents were (2.498, 0.500, 1.000). 
These polarization orbitals were added to the 
Dunning-Huzinaga (10s6p) [20,41,42] basis with the 
HI s orbitals scales to 1.2 which produced a lower 
energy than that of a 6-3 1G** basis using the 
GAM ESS program. The bond length of OH + was 
optimized using GAMERS at the 6-3 1G** basis level 
and fund to be 1.0062 A. The same bond length was 
used for OH" since the anion calculation using the 
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usual HFR-SCF process was not feasible and in any 
case the bond length is only slightly longer than that 
in water. The horizon cutoff value of O.O0038 A opti- 
mized for H 2 0 was also used for OH + and OH". The 
spike in the numerical integral routine was optimized 
by fitting the r c cutoff value so as to obtain as near as 
possible the non-relativistic energy of the HF molecule 
as determined from quantum Monte Carlo calculations 
[21]. The dipole moments for the ions are not very use- 
ful since ion dipoles are origin dependent, but they 
were calculated using the center-of-mass as the origin. 

As we see in Table 1, the difference in energy 
between OH" and OH 4 is 0.497621 hartrees (13.54 
electron volts) according to the correlated-SCF calcu- 
lations. It is clear from the standard SCF energy value 
for H 2 0 that this basis is very good, but not quite at 
the Hartree-Fock limit of energy. In addition, the fit- 
ting of the numerical integration spike so as to nearly 
reproduce the total energy of HF is not exact. These 
two artifacts introduce an energy uncertainty of about 
0.0115 hartrees, but this is less uncertain than that of 
the quantum Monte Carlo energy of Luchow and 
Anderson [22]. Note that the iso-dipoles for H 2 0 and 
HF are very close to the experimental values which in- 
dicate that the calculated wavefunctions are of high 
quality as well as the energy values. 

Since the ionization energy of a neutral H-atom is 
13.60 eV and the energy difference of 13.54 eV would 
convert OH" to OH + , a threshold of about 13.7 eV 
should maintain H + in a solution as well as transfer 
(-2e) through an aqueous solution to or from the 
OH"/OH + system. These calculations indicate that 
there may be an enhancement of current flow above a 
potential above 13.7 V across an aqueous cell and that 
the enhanced concentration of HOOH may be measur- 
able above a potential of 13.7 V. It is worth repeating 
that this estimate is possible largely due to the easy 
convergence of the correlated-SCF process for a nega- 
tive ion species; a process which is formally not defined 
under the usual Hartree-Fock-Roothaan process 
[14,38] and most quantum chemists are familiar with 
the difficulty in treating negative ions using the stan- 
dard Hartree-Fock-Roothaan method. 

Admitting that the correlated-SCF equations are a 
parametrized approximation of the Santilli derivation 
of the Hulten potential [6,31] for a bound electron- 
pair, the method has the advantage of easy incorpor- 
ation into an existing Hartree-Fock-Roothaan 
Gaussian basis program merely by subtracting a small 
'correlation integral' from the usual two-electron inte- 
grals. With some thought one should realize that fitting 
the single parameter (Gaussian screening exponent, a) 
to experimental energies and/or quantum Monte Carlo 
results will incorporate another attraction in the form 
of a magnetic dipole interaction between the spin 
moment of paired electrons. Including the magnetic 
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dipole interaction and substituting a Gaussian form 
for the Hulten exponential potential leaves only a 
simulation of the bound electron-pair 'isoelectronium* 
predicted by Santilli. Thus, these results are for a 
model in which the usual HFT-SCF method is cor- 
rected for at least two attractive interactions of elec- 
trons causing them to approach each other as if they 
were 'sticky*; hence the term 'sticky-electron-pair 
model'. 



5. Conclusions 

In the preceding paper [3,30] of this series we have 
presented a covering of quantum chemistry under the 
name of hadronic chemistry and applied the new disci- 
pline to the construction of a new model of molecular 
structures based on the bonding of a pair of valence 
electrons from different atoms into a singlet quasi-par- 
ticle state called isoelectronium. In the preceding paper 
[3] we then applied the model to the structure of the 
hydrogen molecule by achieving results such as: a rep- 
resentation of the binding energy and other features of 
the hydrogen molecule accurate to the seventh digit; 
an explanation of the reason why the hydrogen mol- 
ecule has only two hydrogen atoms; a reduction of 
computer usage at least 1000-fold; and other advances. 

In this paper we have applied the isochemical model 
of molecular bonds to the water and other molecules 
with similar results. In fact, the isochemical model of 
the water and other molecules presented in this paper 
is supported by the following conceptual, theoretical 
and experimental evidence: 

1. It introduces a new strong binding force (which is 
absent in current models) capable to explain the 
strength and stability of the molecules. 

2. It explains the reason why the water molecule has 
only two hydrogen atoms. 

3. It permits a representation of the binding energy of 
the water and other molecules which is accurate to 
several digits. 

4. It represents electric and dipole moments and other 
features of the water and other molecules also accu- 
rate to several digits. 

5. It permits a reduction of computer usages in calcu- 
lations at least 1000-fold; as well as permitting other 
achievements similar to those obtained for the 
hydrogen molecule. 

Moreover, as it happened for the hydrogen molecule 
[3,30], the value of the radius of the isoelectronium 
computed via dynamical equations in Section 3 has 
been fully confirmed by independent calculations for 
the water and other molecules conducted via the 
Gaussian-lobe basis set. 
The emission of electron pairs in superconductivity 
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has been emphasized in the preceding papers. We indi- 
cate here preliminary yet direct experimental verifica- 
tions of the isochemical model of molecular bonds 
offered by the ongoing experiments on photoproduc- 
tion of the valence electrons in helium indicating that 
electrons are emitted in pairs [23]. The systematic rep- 
etition of these experiments specifically for water is rec- 
ommended here. The statistical percentages of electron 
pairs over the total number of emitted electrons would 
then establish whether the isoelectronium is fully or 
only partially stable. 

We should also recall that Shillady's computer pro- 
gram of isochemistry, called SASLOBE, has been 
shown in this paper to be capable of producing an 
essentially exact representation of experimental data 
on the water and other molecules 1000 times faster 
than conventional programs. 

We should finally note that the representation of the 
binding energy, electric and magnetic moments and 
other characteristics of the water and other molecules 
exact to several digits confirms a main result of the 
preceding paper [3,30], namely, the identification of ex- 
perimental evidence on the insufficiency of quantum 
mechanics and the validity of the covering hadronic 
mechanics for the representation of nonlinear, nonlocal 
and nonpotential-nonunitary effects due to deep over- 
lappings of the 'extended wavepackets' of electrons 
with a 'point-like charge structure'. 

These results evidently provide sufficient credibility 
to warrant systematic theoretical and experimental stu- 
dies on the new isochemical model of water and other 
molecules. 

In closing, we would like to indicate, for future 
study, that the new isochemical model of the water 
molecule implies new predictions of theoretical and 
industrial significance. In particular, in a future paper 
we plan to present theoretical and experimental evi- 
dence on an apparent enhancement of the energy con- 
tent of hydrogen when subjected to certain magnetic 
and other polarizations. 

As additional possibilities we note that, in this work 
the correlated-SCF method is used to easily obtain an 
energy for the OH" anion in water, while the OH + 
ion is easily treated in either the standard or modified 
method. The difference in energy between the eight- 
electron OH + system and the 10-electron OH" system 
is found to be 13.54 eV. This represents the energy 
needed to remove (— 2e) from OH~. This indicates 
there may be a threshold for current flow in terms of 
(-2e) as a quasi-particle in aqueous media at 13.6 eV. 
This voltage will also maintain H + in a solution to 
some extent. Organic alkenes in a solution should 
undergo epoxidation followed by solvolysis to diols 
under the conditions of abundant (-2e). 

Another interesting result is that the natural trace 
amounts of HOOH in water may be increased in water 



by merely placing the sample in an intense magnetic 
field. Positive and negative ions will traverse short arc 
segment paths driven by simple thermal Brownian 
motion in a way which will lead to an increase in col- 
lisions of oppositely charged ions. In particular, OH" 
and OH + may undergo collisions more frequently 
leading to an increase in HOOH. 

This additional HOOH should then be available to 
react with alkenes to form epoxides which will then 
hydrolyze in water to form diols. Such diols would be 
much more soluble in water than the original alkenes. 
This leads to the important possibility that merely 
exposing water-insoluble alkenes to water in a mag- 
netic field will lead to a chemical reaction of the 
alkenes to form modified compounds which are more 
soluble in water. In other words, organic oils contain- 
ing some double bonds may be made somewhat more 
soluble in water just by mechanical emulsification of 
the oils in water in an environment of a high magnetic 
field. 

Thus, mixtures of oils and water could be mechani- 
cally agitated in a magnetic field of several Tesla to 
produce new oils which are chemically similar to the 
original oils (assuming a large organic structure) but 
more soluble in water after exposure to the magnetic 
field. There are numerous important applications of 
such a process to *make oils water soluble 1 . 

The industrial applications of the above novel pre- 
dictions are under development at various corporations 
and they will be treated in detail in future papers. 
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